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REMARKS/ARGUMENTS 

Claims 18-30 are currently pending in the above-identified application. No 
claims have been amended in the present response. Applicants respectfully request 
reconsideration of the pending claims in light of the remarks below. 

Rejections under 35 U.S.C. § 102 : 

Claims 1 8, 19 and 22 remain rejected under 35 U.S.C. § 102(b) as being 
anticipated by U.S. Patent 5,338,686 (Hellerstein, "the '686 patent"). The Examiner asserts that 
Hellerstein teaches a method for determining the rate of degradation of an isotopically labeled 
biopolymer comprising adding a stable isotope-labeled subunit (monomer) to a biopolymer pool, 
collecting first and second samples and measuring the relative abundance of monoisotopic and 
isotopomeric peaks using mass spectrometry, calculating the difference between the peaks of the 
first and second samples and determining the rate of biopolymer decay. The Examiner further 
alleges that serial time points were collected and individual mass isotopomers were plotted over 
time to determine the rate of decay. The Examiner concludes firom these alleged actions that 
Hellerstein teaches every element of the claimed invention. 

In addition, the Examiner questions Applicants assertion that the calculation of 
molar excesses is not the same as the determination of abundance. The Examiner alleges that 
Hellerstein calculates relative abundance using mass spectrometiy and further remarks that she is 
unclear what the difference is between molar excess and relative abundance, since both are 
looking at abundance of monoisotopic and isotopomeric peaks. Further, the Examiner returns to 
the Hellerstein's statement that mass spectral peak heights may be expressed as ratios for 
piarposes of determining relative frequencies, but also states that any calculation means which' 
provides relative values for abundance of such isotopomers in a sample may be used in 
describing the data. 
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Applicants must again respectfully disagree with the Examiner rejection. The 
Examiner has questioned the difference between molar excess and relative abundance and the 
relevance of the difference to any patentable differences between Hellerstein and the present 
invention. Molar excess as described and used by Hellerstein describes a calculation of the 
observed distribution of a mass isotopomer in a polymer containing at least two identical 
subunits. In calculating molar excess Hellerstein first isolated each biopolymer containing 1, 2, 
3,4,... etc, (Mo, Ml, M2, . . . etc.) isotopically labeled subimits. Each of the isolated calculated 
values is for a particular isotopomer species. The determined value for each isotopomer species 
is then compared to a probability matrix that is created on assumptions of precursor pool 
percentages, which in the method of Hellerstein resuhs in the calculation of the molar excess for 
each isotopomer species. Once the molar excess is calculated, the frequency of replacement for 
each isotopomeric species is compared at two different time points to calculate the decay of a 
biopolymer. In a 1999 review published by the inventor, (Hellerstein and Neese, Am. J. Physiol. 
276:E1 146-1 170, 1999), state "that subpopulations of molecules must be distinguishable and 
quantifiable" (page El 147, left column, fourth full paragraph, lines 1 and 2). In addition, they 
state that "there must be combinations possible in the molecule analyzed" (page El 147, left 
column, third full paragraph, lines 1 and 2) and that "the polymers must be intact, or with at least 
two subunits present, because the distribution of the isotopomeric species carries the essential 
information." (page El 147, left column, third full paragraph, lines 7 - 10). This characteristic of 
the invention of Hellerstein is recited in column 6, lines 60 through 64. 

In order to try to further clarify the method of Hellerstein, Applicants provide the 
following analogy also used by Hellerstein in publications from his laboratory. The analogy uses 
a necklace containing colored beads to represent a biopolymer comprising a series of subunits. 
In the analogy, specific beads (within a pattern of beads of many colors that make up the whole 
necklace) can either be colored white (non-labeled) or colored black (isotopically-labeled). 
Naturally occurring necklaces (biopolymers) will contain white beads and some percentage of 
black beads where white beads should be found because black beads, occur naturally (natiu-ally- 
occurring stable isotopes found in nature). When the bead precursor pool is enriched with black 
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beads over white beads, there is an increased probability that black beads will be incorporated 
into the necklaces at the location where white beads are usually found. The method of 
Hellerstein looks at individual necklaces, isotopomeric species, using, in the specific examples 
gas chromatography, measuring the number of necklaces containing 1, 2, 3 ... etc. black beads, 
and comparing these data to probability matrices calculated using different proportions of white 
and black beads available for stringing on necklaces (the precursor pool). Further, Hellerstein 
uses the values measured for the mass isotopomeric peaks to determine whether the biopolymers 
are substituted with 1,2,3, etc. isotopically-labeled subunits. 

Fig. IB of the '686 patent provides a description of the Hellerstein method for 
measuring isotopic decay. In the first box, the molar excess for each necklace, isotopomeric 
species, the polymer having one substitution Mi through the polymer having four substitutions 
M4 is calculated at each time point. In the second step the ratio (R) of the calculated molar 
excess to the sum of the molar excesses for each necklace is calculated. Time points and masses 
for which "R" is stable are used to calculate the decay constant (k). The statement of Hellerstein 
does not disclose or suggest the method of the present invention. 

In the method of the present invention, each isotopomer species is not analyzed 
separately, each necklace, but instead the pooled distribution of isotopomers, or necklaces or 
necklace fragments, is analyzed. There is no need to distinguish subpopulations of polymers as 
data is collected on a pooled population containing both labeled and imlabeled polymer (relative 
abundance). This simpler method can be used even when only one subunit is present within the 
polymer being analyzed. In one aspect of the method of the present invention a necklace 
(biopolymer) with a known pattern of white beads is broken down into fragments. Some of the 
fragments will contain white beads that can be substituted with black beads (A fragments), and 
other fragments will not contain any white beads (B fragments). One or more of the necklace 
fragments must contain white beads (or black beads) for our method to work, but unlike 
Hellerstein there is no need for two repeats of an identical subunit to be present in the 
biopolymer. In the method of the present invention, a pool of specific A fragments is analyzed, 
and the mass isotopomeric distribution of the pooled population is measured by mass 
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spectrometry. In order to calculate biopolymer formation or degradation, the mass isotopomeric 
distribution of A fragments that have been formed in the presence of large excesses of black 
beads is compared to the mass isotopomeric distribution of pooled control A fragments 
(containing a mixture of white beads and naturally-occurring black beads). The percentage of 
newly formed necklaces is determined by mathematical comparisons of mass isotopomeric 
abundances of pooled necklace fragments before and after enrichment with black beads. 

As such, Hellerstein does not teach every element of the claimed invention. As 
above, Hellerstein instead teaches a method that calculates the molar excess of each isotopomer 
species and the frequency for that species. This value is compared at different time points. The 
calculation of the molar excess and frequency for each isotopomer does not disclose or suggest 
the present invention. 

The Examiner has also again referred to a comment in Hellerstein that any 
calculation means which provides relative values for abundance of such isotopomers in a sample 
may be used in describing the data as teach elements of the present invention. Applicants 
respectfiilly must again point out to the Examiner that this statement is made in the context of 
calculating the molar excess for each isotopomer species and not the pooled relative abundance 
for all monoisotopic and isotopomeric peaks for an intact or fragment of a biopolymer. The 
present invention provides a simpler method for determining a degradation rate for a biopolymer 
that is not disclosed or suggested by Hellerstein. 

Applicants note that claims 20, 21 and 23 through 30 are objected to as being 
dependent upon a rejected base claim, but would be allowable if rewritten in independent form 
including all of the limitations of the base claim and any intervening claims. 

In view of the above remarks, the Examiner is respectfully requested to reconsider 
and withdraw the rejection of claims 18, 19 and 22 under 35 U.S.C. § 102(b) as being anticipated 
by U.S. Patent 5,338,686 (Hellerstein). 
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CONCLUSION 



In view of the foregoing, Applicants believe all claims now pending in this 



Application are in condition for allowance and an action to that end is respectfully requested. If 
the Examiner believes a telephone conference would expedite prosecution of this application, 
please telephone the undersigned at 206-467-9600. 



TOWNSEND and TOWNSEND and CREW LLP 

Two Embarcadero Center, Eighth Floor 

San Francisco, California 941 11-3834 

Tel: 206-467-9600 

Fax: 415-576-0300 

BWP/jlv 

61228647 vl 



Respectfully submitted, 





Brian W. Poor 
Reg. No. 32,928 
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Hellerstein, Marc K., and Richard A. Neese. Mass 
isotopomer distribution analysis at eight years: theoretical, 
analytic, and experimental considerations. Am. J. Physiol. 
276 {Endocrinol. Metab. 39): E1146-E1170, 1999.— Mass 
isotopomer distribution analysis (MIDA) is a technique for 
measuring the synthesis of biological polymers. First devel- 
oped approximately eight years ago, MIDA has been used for 
measuring the synthesis of lipids, carbohydrates, and pro- 
teins. The technique involves quantifying by mass spectrom- 
etry the relative abundances of molecular species of a poly- 
mer differing only in mass (mass isotopomers), after 
introduction of a stable isotope-labeled precursor. The mass 
isotopomer pattern, or distribution, is analyzed according to a 
combinatorial probability model by comparing measured 
abundances to theor'etical distributions predicted from the 
binomial or multinomial expansion. For combinatorial prob- 
abilities to be applicable, a labeled precursor must therefore 
combine with itself in the form of two or more repeating 
subunits. MIDA allows dilution in the monomeric (precursor) 
and polymeric (product) pools to be determined. Kinetic 
parameters can then be calculated (e.g., replacement rate of 
the polymer, fractional contribution from the endogenous 
biosynthetic pathway, absolute rate of biosynthesis). Several 
issues remain unresolved, however. We consider here the 
impact of various deviations from the simple combinatorial 
probability model of biosynthesis and describe the analytic 
requirements for successful use of MIDA. A formal mathemati- 
cal algorithm is presented for generating tables and equa- 
tions (appendix), on the basis of which effects of various 
confounding factors are simulated. These include variations 
in natural isotope abundances, isotopic disequilibrium in the 
precursor pool, more than one biosynthetic precursor pool, 
incorrect values for number of subunits present, and concur- 
rent measurement of turnover from exogenously labeled 
polymers. We describe a strategy for testing whether isotopic 
inhomogeneity (e.g., an isotopic gradient or separate biosyn- 
thetic sites) is present in the precursor pool by comparing 
higher-mass (multiply labeled) to lower-mass (single- and 
double-labeled) isotopomer patterns. Also, an algebraic correc- 
tion is presented for calculating fractional synthesis when an 
incomplete ion spectrum is monitored, and an approach for 
assessing the sensitivity of biosynthetic parameters to mea- 
surement error is described. The different calculation algo- 
rithms published for MIDA are compared; all share a common 
model, use overlapping solutions to computational problems, 
and generate identical results. Finally, we discuss the major 
practical issue for using MIDA at present: quantitative 
inaccuracy of instruments. The nature and causes of analytic 
inaccuracy, strategies for evaluating instrument perfor- 
mance, and guidelines for optimizing accuracy and reducing 
impact on biosynthetic parameters are suggested. Adherence 
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to certain analj^ic guidelines, particularly attention to concen- 
tration effects on mass isotopomer ratios and maximizing 
enrichments in the isotopomers of interest, reduces error. 
Improving instrument accuracy for quantification of isoto- 
pomer ratios is perhaps the highest priority for this field. In 
conclusion, MIDA remains the "equation for biosynthesis," but 
attention to potentially confounding factors and analytic perfor- 
mance is required for optimal application. 



THE ASSEMBLY AND DISASSEMBLY of polymers Synthesized 
from repeating monomeric units is a central theme in 
biology. Such polymers may be as simple as fatty acids 
synthesized from acetyl-CoA units or as complex as 
proteins synthesized from amino acids or DNA made 
from nucleotides. Other examples include carbohy- 
drates (e.g., glucose from triose units, glycogen from 
glucose, glycoproteins), porphyrins (e.g., chlorophyll, 
heme), and lipids (e.g., cholesterol, triacylglycerols) . 
Biological polymers may be homonuclear (defined as 
containing subunits that are identical), as in fatty 
acids, or heteronuclear (defined as containing more 
than one type of subunit), as in proteins or polynucleo- 
tides. Despite the importance of polymers in the chem- 
istry of living systems, techniques for determining 
their rates of synthesis or breakdown have historically 
been unsatisfactory (1, 9, 18, 19). As a consequence, 
fields as wide ranging as lipid biosynthesis, protein metabo- 
lism, carbohydrate metabolic regulation, and control of 
cell proliferation have been severely constrained. 

In this article, we will provide an update and eight- 
year perspective on a technique that provides a funda- 
mental solution to the problem of measuring polymer- 
ization biosynthesis. Mass isotopomer distribution 
analysis (MIDA) is a technique based on combinatorial 
probabilities and the labeling patterns in intact poly- 
mers that can be said to provide a fundamental "equa- 
tion for biosynthesis." Although MIDA was first pre- 
sented as a systematic approach to polymerization 
biosynthesis only a few years ago (13-15, 20), a number 
of refinements, alternative calculation algorithms, and 
criticisms have been published since then (4, 5, 22, 31, 
32). We will review here the theoretical and practical 
factors that must be taken into account if MIDA and 
related techniques are to be applied successfully. 

BIOLOGICAL BASIS OF MEASUREMENT 
OF POLYMERIZATION BIOSYNTHESIS 
BY ISOTOPE INCORPORATION 

The principle of isotope incorporation techniques for 
measuring polymerization biosynthesis is, on the sur- 
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face, straightforward. In a biological system, polymers 
that are newly synthesized mix into a pool that also 
contains preexisting polymer molecules. The goal of an 
isotope incorporation study is to quantify the fraction of 
molecules in the mixture that were newly synthesized 
during the label incorporation period (i.e., "what's 
new") and the rate at which the total pool of polymers is 
turning over. To determine the newly synthesized frac- 
tion {fj present in the mixture, one must first establish 
exactly how much label is contained in the population 
of newly synthesized polymers. Dilution of this labeled 
population by the population of preexisting, unlabeled 
molecules can then be determined, according to the 
precursor-product relationship (15, 22, 39, 40). 

The major practical difficulty has been establishing 
how much label is contained in the newly synthesized 
population of molecules. There exists no purely physi- 
cal technique for identifying in a mixed population of 
molecules which ones are new and which are not. No 
classical extraction technique can reveal where differ- 
ent molecules in a population came from or how long 
they have been present. The biochemistry of the precur- 
sor-product relationship provides a possible solution, 
however (Fig. 1), because the precursor pool of subunits 
in a cell has a physical reality and can in principle be 
isolated by extraction techniques. 

Serious problems arise when investigators have tried 
to use surrogate monomer pools to represent the isoto- 
pic content of the true precursor pool (p) (9, 19, 34, 37, 
38), however: e.g., plasma amino acids or free intracel- 
lular amino acids to represent the tRNA-amino acid 
precursor pool for protein synthesis, or ketone bodies to 
represent the acetyl-CoA pool for lipogenesis. Complicat- 
ing factors deriving from subcellular or intracellular 
biochemical organization have been shown to affect 
every class of polymer so far examined in detail, 
including proteins (38). lipids (8, 9), carbohydrates (19, 
34), and nucleic acids (18). 

A SOLUTION TO THE PRECURSOR-PRODUCT PROBLEM: 
THE USE OF COMBINATORIAL PROBABILITIES 

MIDA is based on a model of combinatorial probabili- 
ties. Polymerization biosynthesis can be conceptualized 
as a combinatorial process, with monomeric subunits 
from a precursor pool combining into a polymeric 
collection or assemblage. If the monomeric subunits are 
of more than one distinctive type, i.e., labeled and 
unlabeled, then the population of assembled polymers 
will not be of uniform isotopic composition. The poly- 
mers will exist as distinguishable species containing 
varying numbers of the different types of subunits. 
Some species will include no labeled subunits, some 
will include one labeled subunit, some will contain two 
labeled subunits, and so on. The relative proportion of 
each species of polymer is determined by and can be 
calculated from the binomial (or multinomial) expan- 
sion (Fig. lA). The binomial expansion contains two 
variables, the number of subunits in the collection (n) 
and the probability (p) of each subunit being of a 
particular type. Because the number of subunits in a 
biological polymer is constant and known, the sole 



factor determining the relative proportions of each 
polymeric combination (i.e., the quantitative distribu- 
tion of mass isotopomers) is p, the labeling probability 
in the precursor pool. 

The population of intact polymeric assemblages there- 
fore contains information about the precursor pool that 
is not available by analysis of the monomeric units in 
isolation: the combinations of labeled and unlabeled 
subunits in the polymer population, which are mani- 
fested for statistical or mathematical analysis as the 
distribution of mass isotopomers. This is the central 
insight on which MIDA is based. Because each isotopo- 
meric distribution is uniquely determined by p, each 
distribution is characteristic of and capable of revealing 
the unique value of p from which it was assembled. The 
distribution is, moreover, immutable; it is a fingerprint 
that will persist throughout the lifetime of the popula- 
tion, as long as there is no biological discrimination (no 
isotope effect) between species of the polymer and no 
remodeling of the polymer after its original assembly. 

What is the effect of mixing a population of polymers 
assembled from a precursor pool of labeling probability 
p with a population of polymers assembled from an 
unlabeled precursor pool? Mixing of this sort (dilution 
of the polymer pool) is what happens in a biological 
system when a labeling experiment is performed: newly 
synthesized polymers from the labeled pool mix with 
polymers that were present before the experiment 
began. A key mathematical feature of MIDA is that the 
relationships among those polymeric species that con- 
tain labeled subunits (the internal pattern among 
isotopomers) are unchanged by dilution from an unla- 
beled population of polymers (Refs. 13-15; see CENTRAL 
FEATURES OF MIDA SUMMARIZED) . 

CENTRAL FEATURES OF MIDA SUMMARIZED 

The first rule of MIDA is that there must be combina- 
tions possible in the molecule analyzed. At least two 
repeats of a probabilistically identical subunit must be 
present. Metabolic pathways involving other kinds of 
chemical transformations but no polymerization are 
therefore not amenable to the combinatorial approach. 
Polymers studied must also be analyzed intact, or with 
at least two subunits present, because the distribution 
of isotopomeric species carries the essential informa- 
tion. Any maneuver that reduces the population to 
monomeric homogeneity, such as combustion to carbon 
dioxide for isotope ratio measurements or hydrolysis to 
monomeric subunits before analysis, loses the combina- 
torial information and precludes application of MIDA. 

The second rule of MIDA is that subpopulations of 
molecules must be distinguishable and quantifiable. 
Indeed, it is the variations within a population of 
assembled polymers that carry the information crucial 
for MIDA. The notion that there is a homogenous 
precursor pool and a uniform product pool is replaced 
by the notion of subpopulations of precursors (some A, 
some B) and subpopulations of products (of characteris- 
tic isotopomeric composition in quantifiable propor- 
tions). Any analytic modality must therefore be capable 
of discriminating among different polymeric subpopula- 
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Fig. 1 . Mass isotopomer distribution analysis 
(MIDA) principle. A: combinatorial probabili- 
ties determine the mass isotopomer pattern 
in polymers. In this simulation, natural abun- 
dance or 10% labeled pools of a subunit com- 
bine into a polymer of eight subunits. The 
population of each of these pools will contain a 
characteristic distribution of Mo, M. M, and 
so on, molecular species (mass isotopomers). 
These proportions can be represented as a 
frequency histogram of the mass isotopomer 
pattern in the polymer and can be measured 
by mass spectrometry After correction for 
natural abundance, degree of enrichment of 
precursor pool can be calculated by comparing 
measured patterns of mass isotopomer abun- 
dances with those predicted from theoretical 
precursor pool enrichments. B: simple numeri- 
cal example of MIDA principle. C: three- 
dimensional representation of change in frac- 
tional abundance (A/li) of a particular mass 
isotopomer (in this case, M\ of methyl palmi- 
tate) as a function of p and f. A plane is 
extended at values of Mi = 0.05, 0.10, 0.15, 
0.20 and 0.25, demonstrating the family of 
solutions for all combinations of p and /that 
give this value of A/li. Inset a 2-dimensional 
projection of p vs. fin the plane Ay4i = 0.15. 
Note linearity of A/4i vs. f. 
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tions (species) present within the population. This is 
why radioisotopic methods cannot be used: specific 
activity is measured from the total counts and total 
mass of material present, treated as a uniform popula- 
tion; and it is why average mass measurements by 
electrospray ionization-mass spectrometry also cannot be 
used: a "centroid" average mass collapses all of the popula- 
tion variability in the polymer pool into a single value. 

The third essential concept underlying MIDA is that 
dilution of the monomeric (precursor) and polymeric 
(product) pools affects abundance distributions differ- 
ently. Both sources of dilution can alter the relative 
proportion of polymeric species containing no labeled 
subunits vs. labeled subunits, but only dilution in the 
precursor pool can alter the internal quantitative rela- 
tionships among labeled species. It is this differential 
effect on "amount" (proportion of the polymer popula- 
tion containing any labeled subunits) vs. "pattern" (relation- 
ships within the population of labeled polymers) that 
allows independent calculation of p and f, respectively. 

The model just described involves some simplifying 
assumptions, which are discussed in theoretical is- 
sues FOR THE USE OF MIDA. 

In addition to determining p directly, MIDA offers 
several operational advantages over previous isotopic 
techniques for measuring biosynthesis rates (Table 1). 
One analytic consequence of the need for intact combi- 
nations is that sophisticated mass spectrometric tech- 
niques must be used in cases in which biosynthesis of 
high-molecular-weight polymers, such as proteins or 
oligonucleotides, is being measured. 

Definitions 

It is useful to define terms to avoid ambiguity. The 
following definitions will be used here. 

Isotopes. Atoms with the same number of protons and 
hence of the same element but with different numbers 
of neutrons (e.g., H vs. D). 

Exact mass. The mass calculated by summing the 
exact masses of all the isotopes in the formula of a 
molecule (e.g.. 32.04847 for CH3NHD). 

Nominal mass. The integer mass obtained by round- 
ing the exact mass of a molecule. 

Isotopomers. Isotopic isomers or species that have 
identical elemental compositions but are constitution- 
ally and/or stereochemically isomeric because of isoto- 
pic substitution, as for CH3NH2, CH3NHD. and 
CH2DNH2. 

Isotopologues. Isotopic homologues or molecular spe- 
cies that have identical elemental and chemical compo- 
sitions but differ in isotopic content (e.g., CH3NH2 vs. 
CH3NHD in the example above) (36). Isotopologues are 
defined by their isotopic composition; therefore, each 
isotopologue has a unique exact mass but may not have 
a unique structure. An isotopologue is usually com- 
prised of a family of isotopic isomers (isotopomers) that 
differ by the location of the isotopes on the molecule 
(e.g., CH3NHD and CH2DNH2 are the same isotopo- 
logue but are different isotopomers) . 

Mass isotopomer. A family of isotopic isomers that is 
grouped on the basis of nominal mass rather than 



Table 1 . Some advantages and disadvantages of MIDA 



Advantages 

1. No uncertainty regarding accuracy of true precursor pool esti- 
mate (e.g., concerns about compartmentalization of pools): the 
subunits that actually entered the polymer are the subunits 
that determine mass isotopomer abundance distributions. 

2. No need for isolation of putative precursor subunits. 

3. All parameters are derived from a single measurement (mass 
spectrometry of polymeric product). 

4. May be applied on multiple products simultaneously. 

5. Applicable in principle to all classes of polymers that can be 
stable isotopically labeled. 

6. Can be used in decay phase with high-mass isotopomers to 
achieve or confirm true pulse-chase conditions and determine 
turnover of polymer pool. 

7. Provides information about intracellular metabolite flow into 
and through precursor subunits (kinetics of the precursor) . 

8. Can provide information about tissue source of polymers, sub- 
cellular pools, timing or location of biosynthetic events, 
metabolite channeling, or other characteristics of in vivo bio- 
synthetic processes, on the basis of measurements of precursor 
pool enrichments. 

9. Does not typically require major perturbation of the physi- 
ological system studied, such as pool flooding, inhibition of 
synthesis or removal, exogenous labeling and reinjection, and 
the like. 

10. Not affected by isotope discrimination at the level of the pre- 
cursor pool. 

Disadvantages 

1 . May require relatively large precursor pool enrichments for 
some applications, which could alter metabolism due to non- 
tracer effects or be expensive to perform. 

2. Requires analytic accuracy, not just precision, because stan- 
dard curves cannot routinely be generated (dilution in both 
precursor and product pools would need to be simulated). 

3. Requires analysis of intact polymers, or portions containing 2 
or more repeats of a monomeric subunit, which can be analyti- 
cally difficult, especially for large heteronuclear polymers (e.g., 
proteins, polynucleotides). 



isotopic composition. A mass isotopomer may comprise 
molecules of different isotopic compositions, unlike an 
isotopologue (e.g., CH3NHD, '3CH3NH2, CH3'5NH2 are 
part of the same mass isotopomer but are different 
isotopologues). In operational terms, a mass isoto- 
pomer is a family of isotopologues that are not resolved 
by a mass spectrometer. For quadrupole mass spectrom- 
eters, this typically means that mass isotopomers are 
families of isotopologues that share a nominal mass. 
Thus the isotopologues CH3NH2 and CH3NHD differ in 
nominal mass and are distinguished as being different 
mass isotopomers, but the isotopologues CH3NHD, 
CH2DNH2, 13CH3NH2, and CH3I5NH2 are all of the 
same nominal mass and hence are the same mass 
isotopomers. Each mass isotopomer is therefore typi- 
cally composed of more than one isotopologue and has 
more than one exact mass. The distinction between 
isotopologues and mass isotopomers is useful in prac- 
tice, because all individual isotopologues are not re- 
solved using quadrupole mass spectrometers and may 
not be resolved even by using mass spectrometers that 
produce higher mass resolution, so that calculations 
from mass spectrometric data must be performed on 
the abundances of mass isotopomers rather than isoto- 
pologues. The mass isotopomer lowest in mass is repre- 
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sented as Mq, for most organic molecules, this is the 
species containing all '^c, iR, '^O, '^N, and the like. 
Other mass isotopomers are distinguished by their 
mass differences from Mq (Mj, Mz, etc.). For a given 
mass isotopomer, the location or position of isotopes 
within the molecule is not specified and may vary (i.e., 
"positional isotopomers" are not distinguished). 

Mass isotopomer pattern. A histogram of the abun- 
dances of the mass isotopomers of a molecule. Tradition- 
ally, the pattern is presented as percent relative abun- 
dances, where all of the abundances are normalized to 
that of the most abundant mass isotopomer; the most 
abundant isotopomer is said to be 100%. The preferred 
form for applications involving probability analysis, 
such as MIDA, however, is proportion or fractional 
abundance, where the fraction that each species contrib- 
utes to the total abundance is used (see ). The term 
isotope pattern is sometimes used in place of mass 
isotopomer pattern, although technically the former 
term applies only to the abundance pattern of isotopes 
in an element. 

Monoisotopic mass. The exact mass of the molecular 
species that contains all 'H, '^c, '^N, '^O, ^^S, and the 
like. For isotopologues composed of C, H, N, 0, P, S, F, 
CI, Br, and I, the isotopic composition of the isotopo- 
logue with the lowest mass is unique and unambigu- 
ous, because the most abundant isotopes of these 
elements are also the lowest in mass (23). The monoiso- 
topic mass is abbreviated as mo, and the masses of 
other mass isotopomers are identified by their mass 
differences from Wq (mi, m2, etc.). 

Fractional abundances. The abundances of indi- 
vidual isotopes (for elements) or mass isotopomers (for 
molecules) given as the fraction of the total abundance 
represented by that particular isotope or mass isoto- 
pomer. This is distinguished from relative abundance, 
wherein the most abundant species is given the value 
100 and all other species are normalized relative to 100 
and expressed as percent relative abundance. For a 
mass isotopomer A/, 



fractional abundance of = 



abundance 



2 abundance Mj 



where 0 to /3 is the range of nominal masses relative to 
the lowest-mass (Mq) mass isotopomer in which abun- 
dances occur. 

Afractional abundance (enrichment or depletion) = 



I abundance 



2 abundance Mj 



I abundance Mx 



abundance Mj 



where subscript e refers to enriched and subscript b 
refers to baseline or natural abundance. 

Isotopically perturbed. The state of an element or 
molecule that results from the explicit incorporation of 
an element or molecule with a distribution of isotopes 
that differs from the distribution found in nature 



(Table 2), whether a naturally less abundant isotope is 
present in excess (enriched) or in deficit (depleted) . 

Monomer A chemical unit that combines during the 
synthesis of a polymer and that is present two or more 
times in the polymer. 

Polymer. A molecule synthesized from and containing 
two or more repeats of a monomer. 

CALCULATIONS 

A tradition in mass spectrometric applications has 
often been to express quantitative results as relative 
abundances (each species normalized to the most abun- 
dant species, which is given the value 100). In contrast, 
fractional abundances or analogous expressions are 
generally preferable for MIDA, because the method is 
based on combinatorial probabilities, and probabilistic 
events are most directly represented as fractions of the 
total universe of choices possible. 

THEORETICAL ISSUES FOR USE OF MIDA 

Every model is based on assumptions, which may or 
may not describe real biological systems accurately. It 
is useful to consider and ultimately to be able to 
evaluate or correct for potential deviations from the 
simple MIDA model. For all of the simulations per- 
formed here, calculations were carried out by use of the 
computer algorithms described in the appendix. 

Effect of Variations in Natural Abundance Values 
of the Isotopes of Elements 

The contribution to mass isotopomer distributions in 
a polymer from natural abundance isotopes of the 
elements has to be subtracted, or otherwise taken into 
account, for labeled subunits to be quantified (13, 15, 
20, 22). One of the most obvious questions that is asked 
is whether the theoretical natural abundance values 
selected, especially for '^C, could have significant ef- 
fects on the calculations. We have done calculations for 
several polymers while varying natural '^C fractional 
abundance between 1.08 and 1.11%, the range that 
might be present in biological carbon in mammals (6, 
33), in the calculation algorithm. The results are shown 
for palmitate-methyl ester (Fig. 2) . There is very little 
effect on calculated values of /when p > 0.03. The same 
is true for glucose, cholesterol (15), and other molecules 
(not shown). Thus variations in natural abundance 
values of '^C do not have an important effect on 
calculated parameters. 

Effect of Isotope Discrimination 

Isotope discrimination or fractionation at the level of 
the precursor pool is not a problem within a biosynthe- 
tic model based on analysis of combinatorial probabili- 
ties (MIDA) . The MIDA calculation reveals the isotope 
content of the subunits that actually entered a polymer, 
regardless of their relation to the isotope content of 
biochemical intermediates leading to these subunits. If, 
for example, there were a 10% discrimination against 
[^HsJleucine by leucyl-tRNA synthase or a 1 0% discrimi- 
nation against [^Hs] leucyl-tRNA by the ribosomal pro- 
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Fig. 2. Effect of varying the natural abundance '^C 
value on mass isotopomer abundances and MIDA calcu- 
lations for methyl palmitate. Values of /are calculated 
using our standard algorithm ('^C = 1.09%), when 
sample in fact has natural abundance '^C = 1.08, 1.10, 
and 1.11%. Calculated values of f are shown as a 
function of p (calculated iteratively, at intervals of 
0.001). Inset, blown-up region of graph between p = 0.05 
and 0.15. 



0 0.05 0.10 0.15 0.20 0.25 

P 



tein synthesis machinery, it would not affect the calcu- 
lations of for true p; it is the ^H^ enrichment of the 
leucine subunits that actually entered the protein that 
determines the mass isotopomer abundance pattern. 
This pattern will reveal the true PH3] leucine precursor 
subunit enrichment for biosynthesis, even if this value 
is different from the tRNA-leucine or free leucine 
enrichments, and the calculated fractional synthesis 
contribution would be correct. Measurements of p, f, 
and the like would then be accurate with MIDA but not 




0 10 20 30 40 50 60 70 80 90 100 



True f (%) 

Fig. 3. Loss of linearity in the relationship between measured A/l,, 
and true f, especially at high values of p, if an incomplete ion 
spectrum is sampled. In this example, only ions Mo to Mz are 
measured in palmitate-methyl-ester at various values of p. The ratio 
A/liM7 (measured excess M/asymptotic M values) should reflect f 
and rise linearly as /increases from 0 to 100% (e.g., see Fig. 1 Q. As p 
increases and the proportion of unmonitored ions rises, however, 
linearity is lost and true / is underestimated. Simulations were 
shown at data points of / 



if tRNA-leucine were used, in this example. Thus 
combinatorial probability analysis is unique among 
isotope kinetic approaches in that its validity is not 
altered by isotope discrimination during biosynthesis. 
In contrast, isotope effects on the metabolism of the 
polymer once synthesized (e.g., effects on clearance) 
will affect kinetic measurements, because the behavior 
of labeled polymers will not reflect that of the general 
pool. 

Effect of Incorrect Value for Expected Number 
ofMonomeric Subunits 

Lee et al. (22) demonstrated elegantly that n (the 
number of precursor subunits actually present in a 
polymer) can be determined experimentally by using 
the same principles of probability analysis that are 
used for determining / Instead of a reference table for p 
vs. mass isotopomer pattern at a known value of n, one 
can generate a reference table for n vs. mass isotopomer 
pattern at a known value of p. The true value of n can 
then be inferred from the experimental data. This 
technique is possible only when there exists an indepen- 
dent method for determining p; the measurement of 
body water enrichments during ^HaO incorporation 
experiments represents a unique situation that permit- 
ted this application (22). 

Calculation of p and f When the Complete 
Ion Spectrum Is not Sampled 

For any number of reasons, the mass spectrometrist 
may choose not to monitor all of the ions in a mass 
isotopomer envelope (e.g., for convenience, to maximize 
dwell time on the most abundant ions, or to avoid 
contaminating ions) . An important property of combina- 
torial probabilities is that the calculation of p is not 
affected by the choice of ions selected for monitoring. 
The internal pattern and relationship among excess 
mass isotopomers are fixed and characteristic of p and 
n regardless of which particular masses are monitored. 
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As long as the appropriate equation is used for the 
masses under consideration, the choice of masses moni- 
tored will not influence calculation of p. Surprisingly, 
this is not the case for calculation of f, which is affected 
by incomplete ion spectrum sampling. This is because 
of a somewhat unexpected mathematical feature of 
mixtures of numerical distributions (e.g., populations 
of mass isotopomers) : dilution is not linear when the 
proportion of the total population monitored is different 
in the natural abundance and enriched populations. As 
noted in the APPENDIX (and see Fig. 1 C), the mathemati- 
cal object of solving for p in step-wise calculation 
algorithms is to linearize the relationship between 
abundance of a particular mass isotopomer (A^ and the 
molar fraction of its associated molecular population in 
the mixture [f), so that / can be solved algebraically 
from Ax. When different proportions of the total ion 
envelope are monitored for different populations, for 
example, as occurs when a high p generates high mass 
isotopomers that are not monitored, the linear relation- 
ship between any Ax and / (fractions of the mass 
isotopomer and molecule in the population, respec- 
tively) is lost (Fig. 3). Stated in intuitive terms, when 
higher masses are not monitored, a mole of isotopically 
enriched molecules will contribute fewer ions to the 
total spectrum sampled than a mole of natural abun- 
dance molecules. The molecular mixture is thereby 
weighted in favor of the more completely sampled 
molecular population (the unlabeled population), and a 
correction has to be made to put equal weight on each 
molecular population in the mixture. 

A numerical example follows. If theoretically 
represents 20% of the Mq through M2 ions in natural 
abundance molecules and 40% of the Mq through Mz 
ions at p = 0.10, but only 90% of the envelope is 
contained in Mq, through M2 at p = 0.10 compared with 
100% in Mo through Mz for the natural abundance 
molecules, what would be the effect of mixing these 
populations and monitoring only Mq through A/2? It 
should be apparent that the enriched population will 
not contribute the theoretical 40 ions of Mi for every 20 
ions of Ml from unlabeled molecules in an equimolar 
mixture when only Mq through M2 are measured, but 
will contribute only 36 (= 40% X 90% of the mole in the 
envelope monitored) for every 20 (= 20% X 100% of the 
mole of natural abundance isotopomers). It would 
therefore be a mistake to use percentages of the ions 
monitored to represent percentages of the entire popu- 
lation when mixing populations with different percent- 
ages of ions monitored, because /will be systematically 
underestimated. The solutions to this confounding fac- 
tor are straightforward: either monitor an essentially 
complete ion spectrum for the molecules under consid- 
eration or include a mathematical correction for un- 
equal ion spectrum sampling in the calculation algo- 
rithm. 

An algebraic correction is derived and presented in 
the APPENDIX (Eq. A 9b) for instances of significantly 
incomplete ion spectrum monitoring. This equation 
corrects for the proportion of ions monitored at the 
measured value of p present relative to the proportion 



monitored in unlabeled molecules. By use of this correc- 
tion factor, mixtures of labeled and unlabeled molecules 
are again reduced to linear combinations of mass 
isotopomers, from which dilution of molecules can be 
calculated simply. This correction is generally ex- 
tremely small and has no practical impact on most 
calculations (Fig. 4), because >98% of the ions within 
an isotopomeric envelope are typically monitored for 
most labeled molecules. Failure to consider the effects 
of incomplete ion spectrum sampling can contribute to 
underestimation of values of fin special cases, however, 
such as very high values of p, if high masses are not 
monitored. The impact of incomplete ion spectrum 
monitoring has not to our knowledge been considered 
previously or corrected for in other MIDA calculation 
algorithms. 

Effect of Analytic Artifacts — Fragment Ions 

An important analytic concern is the effect on calcu- 
lated parameters of fragment ions generated in the 
mass spectrometer source. The most common problem 
is the M-i fragment, which can come from extraction of 
H from the parent ion. The effect of contaminating M_i 
fragments can be simulated using the calculation algo- 
rithm described here (see appendix). Contaminating 
M_i fragments representing 0, 1, and 5% of ions were 
included in measurements of palmitate synthesis from 
acetyl-CoA (not shown). The contaminated fractional 
abundance distribution of mass isotopomers was calcu- 
lated by overlaying the weighted intact and fragment 
ion distributions, and calculations of p and f were 
derived as though they had been generated as experi- 
mental data. Even at 5% contaminant contribution, 
calculated /"was >99%. Because the fragment contrib- 
utes to natural abundance distributions as well as 
enriched distributions, its impact is substantially cor- 
rected through the subtraction of natural abundance 
values. It is nevertheless worthwhile to look for contami- 
nating fragment ions in the mass spectrum. 




Fig, 4, Error in calculation of / resulting from incomplete ion 
spectrum monitoring of different degrees. Three different values for 
true f{25. 50, 75%) are shown. Note that calculated /remains a95% 
of actual value ais long as —85% of the ion spectrum is monitored. 
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Effect oflsotopic Disequilibrium in the Precursor Pool 

Evaluating the theoretical impact of isotopic disequi- 
librium. A biosynthetic system may or may not exhibit 
isotopic equilibrium among different pools of a mono- 
meric subunit from which a polymer is synthesized. 
Accordingly, there may not always be a single value for 
p. A physiologically relevant example of a biosynthetic 
system that does not necessarily have a single value for 
pis gluconeogenesis. There are actually two precursors 
comprising the gluconeogenic triose-phosphate pool, 
dihydroxyacetone phosphate (DHAP) and glyceralde- 
hyde-3-phosphate. Isotopic equilibrium between the 
triose-phosphates may not always be complete. One 
can simulate the effect of various degrees of isotopic 
disequilibrium between DHAP and glyceraldehyde-3- 
phosphate and determine the extent to which calcu- 
lated values of p and / would be distorted if the 
standard MIDA reference table were applied to the 
Afractional abundances generated (27). If we allow the 
average pto range between 0.05 and 0.15 and vary pin 
pool 2 from 1.0 to 2.0 times the pool 1 value, the 
consequences for the /can be calculated. When p in pool 
1 = p in pool 2 (i.e., when isotopic equilibrium is 
present) , the calculated value of /is of course exactly as 
expected from the MIDA tables (100%), whatever the 
value of p. When p in pool 2 is double that of pool i, 
calculated /is —12% higher than the true value for all 
values of p; when pool 2 is made 50% above pool 1, /is 
within 4% of the true value. An interesting general 
mathematical result to emerge from this analysis is 
that / is always overestimated (>100% the actual 
value) if isotopic disequilibrium exists within precursor 
subunit pools (27). The practical implication for measur- 
ing gluconeogenesis in particular is that isotopic disequi- 
librium in the triose-phosphate precursor pool is in 
principle unlikely to represent a major problem, be- 
cause the MIDA calculations will work well unless 
DHAP and glyceraldehyde-3-phosphate are differen- 
tially enriched by a factor of >2 (27). Similar calcula- 
tions can be applied to other polymers of interest. 

Correcting for documented isotopic disequilibrium. If 
one is not satisfied with theoretical arguments discount- 
ing the importance of isotopic disequilibrium within a 
precursor pool, it is possible to modify the calculation 
algorithm to incorporate deviations from isotopic equi- 
librium within precursor pools. Theoretical tables can 
be generated by modifying the algorithm described 
above, if the degree of isotopic disequilibrium is measur- 
able. In the case of gluconeogenesis, for example, mass 
spectrometric fragmentation of the molecule into "top" 
(C-1 to C-3) and "bottom" (C-4 to C-6) halves can reveal 
whether labeling was equal in DHAP and glyceralde- 
hyde-3-phosphate, respectively (26). If enrichments 
differ by a certain proportion, then this value can be 
incorporated into the probability calculations to gener- 
ate an appropriate, individualized standard curve that 
adjusts asymptotic values appropriately. Thus individu- 
alized standard curves can be generated for each 
experiment on the basis of the observed degree of 
isotopic disequilibrium within precursor pools, if the 



latter is significant and can be determined experimen- 
tally. Because the mathematical approach that we have 
described (see appendix and Refs. 13-15) is based on 
empirical relationships between derived values (Afrac- 
tional abundances) rather than expressions of the pure 
binomial or multinomial expansion, deviations from 
the simple combinatorial probability model can be 
accounted for relatively easily and without compromis- 
ing mathematical rigor. 

Effect offsotopic Inhomogeneity in the Precursor Pool 
(e.g.. More Than One Biosynthetic Site, an Isotopic 
Gradient Across a Tissue, or Time Variations in p) 

A potentially more important deviation from the 
simple MIDA model is if the polymer is made in more 
than one anatomic location and p is not equal in each 
site, or if the value of p changes over time. This 
situation differs from the situation of isotopic equilib- 
rium within a precursor subunit pool (see Effect of 
Isotopic Disequilibrium in the Precursor Pool) . Each of 
the polymer populations synthesized has a single pre- 
cursor enrichment, but there is more than one pool of 
polymers present, whereas in the former case, each 
polymer molecule has more than one precursor enrich- 
ment but there is only a single pool of polymers. 
Examples of more than one anatomic location for 
biosynthesis might include extrahepatic and hepatic 
gluconeogenesis or cholesterogenesis, or labeling gradi- 
ents within a precursor pool across a tissue. The same 
considerations would apply if p changed over time 
during an experiment; polymers made at different 
times would have dilTerent mass isotopomer patterns. 

The consequence of these scenarios is that, instead of 
a single binomial distribution, there will be a mixture of 
distributions from the different values of p in newly 
synthesized polymers. This mixture of distributions 
itself mixes with, and needs to be distinguished from, 
the natural abundance distribution that represents old 
molecules. Any attempt to model biosynthesis as two 
populations or two distributions, enriched and unen- 
riched, will not be rigorously correct mathematically 
but will instead be an approximation, because binomial 
or multinomial expansions cannot themselves be aver- 
aged (i.e., are not linear) (21). The Mz isotopomer 
changes approximately as the power of 2 for changes in 
p, the M3 as the power of 3, etc., so that different 
expansions cannot be combined and averaged as if they 
were linear. 

The practical questions are, how much does this 
matter? what impact on estimated parameters will 
there be if p is inconstant in time or space? and can it be 
identified or corrected for when present? A simulation 
for gluconeogenesis in two tissues or at two time points 
with different p values has been presented elsewhere 
(26): /remains >0.8-0.85, even when pool 2 enrich- 
ment is 2-3 times pool 1. Another example is a gradient 
in precursor pool enrichment across a tissue (Fig. 5 A). 
If one models 10 pools contributing equally to gluconeo- 
genesis with a labeling gradient that spans an approxi- 
mately twofold range [e.g., 0.105-0.195 molar excess 
(ME)], the consequence is minor (underestimation of 
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Fig. 5. Effects of isotopic inhomogene- 
ity in the precursor pool on calculated 
values of p and f. A: effect on gluconeo- 
genesis of an isotopic gradient within a 
tissue. Ten sites are simulated to con- 
tribute equally to gluconeogcnesis, with 
different degrees of isotopic gradient in 
the region around p= 0.15 (from 0.125- 
0.175 up to 0.01-0.29, at the intervals 
stated). B: use of high-mass isotopomer 
analysis to identify inhomogeneity in 
the precursor pool, or in this case, two 
precursor pools. Isotoponieric enve- 
lopes resulting in palmitate-methyl- 
ester if two biosynthetic sites contrib- 
ute equally, one at p = 0.05 and one at 
p = 0.30, are shown (7eA). The envelope 
from the summed distribution is then 
compared with the envelope expected 
forp = 0 . 1 75 (the average va lue) (right) . . 
C. use of high-mass isotopomer analy- 
sis to identify inhomogeneity in the 
precursor pool, in this case, an isotopic 
gradient in the precursor pool. A label- 
ing gradient from 0.03 to 0.30 for p is 
compared with a homogenous pool at 
the value of p (0.166) calculated from 
the observed Mi-to-Mi ratio. 
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true values of /by a factor of <5%, Fig. 5A). Even for a 
fourfold gradient (0.06-0.24 ME), fis only underesti- 
mated by -15%; i.e., if the actual value of / = 0.50, 
measured /will be 0.425. Only at very large gradients 
(e.g., 15-fold, from 0.02 to 0.29 ME) is even a 25% 
underestimation observed. An analogous situation can 
be simulated for lipogenesis, with an isotopic gradient 
of acetyl-CoA across a lipogenic tissue such as liver. If 
10 pools contribute equally to lipogenesis, with a gradi- 
ent from 0.03 to 0.30 ME (at intervals of 0.03), /is 87.3% 
instead of 100% (i.e., underestimated by 12.7%). If a 
gradient from 0.02 to 0.10 ME is simulated, with 100 
pools contributing equally to lipogenesis, the value of / 
calculated by MIDA by use of Mi and M2 isotopomers is 
90. 1%. Thus inconstancy in p over time or space means 
that the simple binomial model becomes an approxima- 
tion rather than an exact description of biosynthesis, 
but the practical impact varies according to physiologi- 
cal conditions and typically is fairly small. An investiga- 
tor is able to evaluate the likelihood of significant error 
and the practical acceptability of this degree of error by 
performing a simulation of this type. 

It is also possible to identify, and even correct for, 
isotopic gradients or variations in precursor pool enrich- 
ments by application of combinatorial principles. If a 
large isotopic gradient is present or two precursors of 
different enrichment contributed to biosynthesis, there 
will be a divergence between the isotopomer pattern in 
high-mass (multiply labeled) vs. low-mass (single- and 
double-labeled) polymers (Fig. 5, 5 and C). If a gradient 
exists within the precursor pool, high-mass species will 
be produced that would otherwise never be observed if 
the average value of p were uniformly present. The 
pattern of higher-mass isotopomers predicted by analyz- 
ing lower masses will not be observed; similarly, the 
pattern of low-mass isotopomers predicted by the higher 
masses will not be met (Fig. 5, Band C). 

The occurrence of "inappropriate" multiply labeled 
species relative to the pattern among the less-labeled 
species (Fig. 5, Sand Q can therefore be used diagnos- 
tically to confirm or exclude the existence of an isotopic 
gradient. One simple approach is to monitor higher- 
mass Isotopomers and compare calculated values of p 
from higher- vs. lower-mass relationships. In the case 
of a labeling gradient across a lipogenic tissue spanning 
0.03 to 0.30 (Fig. 5C), for example, the pattern of excess 
M3/M2 isotopomers in palmitate indicates p = 0.213; 
the ratio in M4/A/3 indicates p = 0.226, whereas Mz/M, 
isotopomers reveal a lower p (0.166). The divergence is 
even greater for simple two-precursor-pool systems 
(Fig. 5B).lf there were an equal mixture of a palmitate 
population derived from 0.30 and a 0.50 value of p, 
analysis of Mi and M2 would reveal p = 0. 173, whereas 
analysis of M3 and M4 would indicate p = 0.29. In 
contrast, if a single, homogenous precursor pool is 
present, calculated values of p are identical whatever 
masses are used for its calculation (15). The finding of 
different p values by analysis of high- vs. low-mass 
isotopomer patterns represents a technique for identify- 
ing the presence of an isotopic gradient. 



Conversion of Fractional Synthesis Values into 
Chemical Fluxes: Combining MIDA Calculations 
with Administration of Exogenous Stable 
Isotope-Labeled Polymers 

Expression of synthesis as rates in chemical units 
(mass/time) requires an estimate of the turnover of the 
polymer pool being sampled in addition to the fraction 
of the polymer pool that came from endogenous synthe- 
sis during the time period studied. Probability consider- 
ations demonstrate that high-mass isotopomers are 
uniquely useful in the labeled-polymer decay phase, 
because problems from persistent isotope incorporation 
are avoided for multiply labeled species even if pulse/ 
chase conditions do not exist (i.e., if precursor subunits 
continue to contain a low level of labeling because of 
isotope recycling or slow turnover of the precursor 
pool). This application of MIDA has been discussed in 
detail previously (15, 25). The turnover of the polymer 
can also be measured by analyzing the rate of rise 
toward plateau during the label incorporation phase, as 
discussed elsewhere (16). Alternatively, exogenously 
labeled polymers can be administered (to measure 
turnover by dilution) concurrently with a biosynthetic 
incorporation experiment. For this last approach, how- 
ever, potential interference by the exogenous label with 
the isotope incorporated via biosynthesis has to be 
accounted for. 

For example, it is useful to measure the plasma 
glucose turnover concurrently with fractional gluconeo- 
genesis to determine the absolute rate of gluconeogen- 
esis. If [l-^HJglucose or [6,6-^H2]glucose is used for 
turnover, labeled species have the same nominal mass 
as the isotopomers analyzed for MIDA calculations (Mi 
and M2). Because quadrupole mass spectrometers do 
not have sufficient resolving power to distinguish be- 
tween and '^C exact masses, a technique for deter- 
mining the contributions from [^HJglucose vs. ['^CJtri- 
ose-phosphate is required. This can be achieved by 
analyzing a derivative that is stripped of the labeled 
hydrogen (e.g., aldonitrile-pentaacetate or saccharic 
acid derivatives from which position 1 and 1,6 hydro- 
gens are removed, respectively) in addition to a deriva- 
tive that contains both inputs (e.g., pentaacetate). The 
MIDA calculation on the derivative stripped of is 
routine. The problem, however, is how to "subtract" or 
correct the '^C contribution from the combined spec- 
trum to establish the labeling, by difference. 

A calculation algorithm can be used to correct for the 
underlying isotopomeric distribution from incorpora- 
tion of gluconeogenic precursors, for example, to 
measure [^H] glucose enrichment (7). The glucose mol- 
ecules present during a simultaneous measurement of 
fractional gluconeogenesis and glucose turnover consist 
of a mixture of three populations: gluconeogenic prod- 
uct molecules arising from the labeled triose-phosphate 
precursor pool, labeled glucose molecules infused exog- 
enously as tracer, and unlabeled molecules with a 
natural abundance distribution. The key is that the 
isotopomeric distribution of each of these components 
is known. The infused tracer ([6,6-^H2]- or [l-^H]glu- 
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cose) has an isotopomer distribution that is easy to 
calculate; natural abundance glucose also has a known 
distribution; and the gluconeogenic population has a 
distribution of isotopomers that is a function of p and 
that is measurable from the deuterium-stripped deriva- 
tive. The distributions from each of these three compo- 
nents of the mixture can therefore simply be added to 
construct a theoretical standard curve, simulating the 
effect of adding ^H-labeled glucose to mixtures of the 
other two populations at the measured f. If gluconeogen- 
esis / = 0.33 and p = 0.15, for example, the calculation 
consists of adding 0.33 times the abundance of each 
isotopomer from gluconeogenesis at p = 0.15, then 
(0.67 - z) times the abundance of each isotopomer from 
natural abundance glucose, and z times the abundance 
of each isotopomer from the [^H] glucose, where z is the 
fraction of ^H-labeled glucose added to generate the 
theoretical standard curve. Two or more values of zare 
simulated to construct a linear standard curve, wherein 
the isotopomers of interest are plotted against z to 
generate a slope and intercept for calculation of 
enrichment and dilution in the intact molecule (Fig. 6). 
This algorithm must be applied separately for each 
time point sampled, because each sample will have a 
unique p and /(8, 36), and thus a unique slope and 
intercept for the standard curve. 

Assessing the Sensitivity of MIDA-Calculated 
Biosynthetic Parameters to Measurement Error 

To plan isotopic experiments and evaluate results, it 
is useful to know how sensitive the derived parameters 
are to analytic imprecision or inaccuracy. Error- 
sensitivity analysis can be performed by using the 
MIDA calculation algorithm presented here. An ex- 
ample with gluconeogenesis has been presented previ- 
ously (26). This analysis reveals that an analytic coeffi- 
cient of variation of 1% in estimates of the ratio of A/I2 
to AAi at p = 0.15 alters estimates of pby ~2% and /by 
1.6%. In contrast, at p = ~0.05, an analytic coefficient 
of variation of 1% alters estimates of p by 5.2% and /by 
4.9% (26). It is apparent from this analysis that the 
experimenter is better advised to aim for p in the range 
of 0.15 rather than 0.05 to reduce sensitivity of final 
results to analytic error. Because analytic precision of 
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Fig. 6. Standard curve for calculation of [6,6-^H2]glucose enrichment 
in the presence of an underlying '^C combinatorial distribution from 
gluconeogenesis. In this example, the value forpis 0.15. 



better than ±0.02 mole percent excess (MPE) can be 
attained by using multiple replicate analyses, several 
time points, a- and (3-anomers of glucose, and other 
analytic strategies (26) , this analysis also reaffirms the 
reproducibility of estimates by MIDA for applications 
such as gluconeogenesis, when pis in the range of 0.15. 
The interaction between enrichment achieved in a 
labeling experiment (i.e., signal) and sensitivity to 
analytic error is discussed further in ANALYTIC AND 
EXPERIMENTAL DESIGN CALCULATIONS. 

ALTERNATIVE CALCULATION ALGORITHMS 

MIDA as we now understand it was developed from 
1990 to 1992 by Hellerstein and co-workers (13-15, 17) 
and Kelleher et al. (20), working independently. MIDA 
is defined by its purpose, its method, and its calcula- 
tions. The purpose or field of MIDA is to measure the 
synthesis of biological polymers, the isotopic labeling of 
the monomeric precursor pool from which the polymers 
were assembled, and related kinetic parameters. The 
method involves measurement and analysis of mass 
isotopomer abundance distributions in intact polymers 
according to a combinatorial probability model, after 
introduction of a stable isotopically labeled monomeric 
subunit. The calculation approach uses the binomial or 
multinomial expansion as a basis for interpreting incor- 
poration of isotopically labeled repeating subunits in 
the intact polymer and for inferring dilution in the 
monomeric and polymeric (precursor and product) pools. 

It is in the area of a calculation algorithm that 
modifications have been presented (1, 5, 22, 32) since 
the original approaches described by Hellerstein and 
co-workers (13-15, 17) and Kelleher et al. (21). All the 
calculation algorithms presented so far, however, share 
fundamental assumptions and postulate a common 
model of biosynthesis. All postulate a combinatorial 
(binomial/multinomial) precursor-product biosynthetic 
model, and all postulate two confounding factors that 
may then modify the simple binomial/multinomial distribu- 
tions: first, the natural abundance isotopes in the molecule, 
which interact with label-derived isotopomers; and second, 
the fact that two sources of dilution exist in biosynthetic 
systems (in the product pool as well as the precursor pool) 
and that these two sources of dilution influence isotopomer 
abundances differently in the product. 

Only two general solutions have been proposed for 
each of the problems just noted. The influence of 
natural abundance isotopes on label distributions is a 
strictly computational problem. The solution has been 
either to create a computational model that incorpo- 
rates all sources of isotope from both labeled precursor 
and natural abundance isotopes, and thus does not 
conform to a simple binomial distribution (Hellerstein 
and Neese, Ref. 15, and Kelleher et al., Ref. 20), or to 
transform the data in a way that removes the influence 
of natural abundance isotopes and restores a pure 
binomial distribution from the labeled precursor (Lee 
et al., Ref. 22, and Chinkes et al., Ref. 5). The problem of 
two sources of dilution, in contrast, reflects a biological 
issue. Again, there have been two computational solu- 
tions proposed. Most methods (Hellerstein, Lee, and 
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Chinkes) have used a step-wise approach. Ratios among 
natural abundance-corrected terms are first computed. 
Use of internal ratios of isotopomers in the polymer 
sample analyzed removes the effect of varying product 
dilution, because all isotopomers sampled in the la- 
beled molecules are equally diluted by natural abun- 
dance molecules, so that precursor pool dilution can be 
calculated independently. Once this unknown (p) is 
solved, the second unknown {/} can be solved algebra- 
ically (Fig. IQ. Alternatively, the two sources of dilu- 
tion (Fig. 1 C) can be solved simultaneously, by best fit 
(nonlinear regression analysis) of multiple solution sets 
for the two unknowns taken together (21). 

The most important point is that all the calculation 
algorithms presented to date give essentially identical 
results. Thus data of Byerley et al. (4), from ^HaO 
incorporation into cellular cholesterol in cultured cells, 
give identical values for p and /whether calculated by 
the method of Lee (22) or by our approach (unpublished 
results). Identical results are obtained when data are 
analyzed by the approach of Kelleher et al. (21) and 
ours (unpublished observations; and T. Masterson, 
personal communication, April 1994). Chinkes et al. (5) 
also compared calculation algorithms using simulated data 
and concluded that results were essentially identical. 

It should be noted that all these approaches are of 
roughly equal computational complexity. They all re- 
quire a computer program', to calculate abundances 
and to generate either a model or algebraic correction 
factors, and a specific software package that has to be 
applied separately for each molecule analyzed. 

ANALYTIC AND EXPERIMENTAL DESIGN 
CONSIDERATIONS 

The single most difficult problem facing the use of 
MIDA at present, both in theory and in practice, relates 
to quantitative accuracy of measurements, i.e., the 
analytic performance of mass spectrometers. This prob- 
lem is widely recognized by workers in the field but has 
only rarely been noted in the literature (3, 10, 11, 23, 
24, 26, 30). MIDA is based on analysis of numerical 
distributions in the context of a model of combinatorial 
probabilities. If the instrument generates inaccurate 
numbers, measured distributions will no longer reflect 
the actual isotopomeric distributions present. The ac- 
tual effect on kinetic estimates (p, f) will depend on the 
nature and extent of the experimental inaccuracy. How 
then does one apply equations that are based on a 
model of combinatorial probabilities if the numbers do 
not fit the actual distributions generated? The best 
solution would be if the experimenter understood the 
cause and exact nature of instrumental inaccuracy: if 
one knew that the explanation was, for example, inad- 
equate resolution of adjacent masses by the mass 
analyzer, a suitable correction algorithm might be 
applied. This is in essence what is done with liquid 
scintillation counting of radioisotopes, wherein and 



' A computer program for calculating theoretical tables and spread- 
sheets is available and wiU be sent by the authors on request. 



'■•C spillover is accounted for, so that each isotope can 
be independently measured. Unfortunately, the ana- 
lytic basis of quantitative inaccuracy by current mass 
spectrometers is not understood in a sufficiently defini- 
tive way at present to allow simple correction (see 
Strategies for Evaluating Quantitative Instrument Per- 
formance and Data Acceptability) . Another reasonable 
approach is to use standard curves, as one does when 
measuring dilution of an exogenous labeled product. 
There are problems also that make this difficult for 
biosynthetic MIDA methods, however. These we will 
discuss. 

Several questions concerning mass spectrometric 
quantitative inaccuracy need to be addressed. We will 
not review these issues extensively here but will note 
some practical implications of each. 

The Nature of Quantitative Inaccuracy in Mass 
Spectrometric Measurement of Isotope Ratios 

Surprisingly little literature exists concerning the 
mass spectrometric causes of deviations between ex- 
pected and measured abundances of mass isotopomers 
in organic analytes (10, 11, 30). The abundance of one 
mass relative to another might be overestimated by a 
fixed proportion, by a fixed amount, or by random error. 
Each type of error would have different implications for 
attempted correction algorithms. One also needs to 
know whether instruments drift over time, so that 
standards analyzed near to a sample in time will reveal 
the existence of error in the sample, or whether error 
occurs in an erratic, unpredictable way. Unfortunately, 
few data and no consensus exist on these questions, 
although several mechanistic possibilities can be consid- 
ered. 

One potential cause of inaccuracy is incomplete 
resolution of adjacent ions in the ion envelope (peak 
tailing), resulting in contamination of adjacent mass 
channels. If the mass analyzer is operating at the limit 
of its mass-resolving capacity, the degree of misidentifi- 
cation of ions due to ion scattering and peak tailing (24) 
could vary from run to run. According to this explana- 
tion, mass analyzers that achieve significantly better 
resolution (e.g., magnetic sector compared with quadru- 
ple mass analyzers) would be predicted to exhibit 
better accuracy. This prediction has not yet been system- 
atically tested. Another prediction is that higher abun- 
dances in adjacent mass isotopomers should worsen 
resolution (increase peak tailing), resulting in worse 
concentration sensitivity of fractional isotopomer abun- 
dances. Abundance sensitivity (24), the observation 
that the most abundant ion in an envelope tends to be 
underestimated quantitatively, has been identified as a 
problem in the field of isotopic analysis of elements by 
use of isotope ratio mass spectrometers. The physical 
explanation and methods for instrumental correction 
continue to be debated, however. A second possible 
cause of inaccuracy is nonlinearity of the detector 
response at different abundances or for different ions. If 
the detector output at each mass does not faithfully and 
consistently reflect the number of ions that reach it, 
numerical distributions will be skewed (10). Detector 
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Table 2 . Instrument factors that might influence measured abundances of mass isotopomers 



Direction of Instrument Dependence 

Dependent 

Effect Effect Effect on Amount 

Cause Effect forM„ forM, for of Sample? TOP Quadrupole Sector 



Bimolecular ion chemistry 


Addition of H 


- 


+ + 


+ 


Yes 


Yes 


Yes 


Yes 


Bimolecular ion chemistry 


Addition of 2 H units 




NA 






Yes 


Yes 


Yes 


Isotope effects during fragmentation 


Misrepresentation of relative abun- 


Var 


Var 


Var 


No 


Yes 


Yes 


Yes 


Incomplete resolution (peak tailing) 


Misidentification of adjacent ions 






-1- 


No 


Yes 




No 


Decreasing ion transmission effi- 


Mass-dependent MS response 


+ 






No 


No 


Yes 


Some 


ciency with mass of ion 


















Nonlinear signal amplification 


Nonlinear signal output 




+ + 




Yes 


Yes 


Yes 


Yes 


Saturation of multiplier 


A4o underrepresented 








Yes 


Yes 


Yes 


Yes 


Velocity dependence of multiplier 


Mass-dependent MS response 


-1- 






No 


Yes 


Yes 


Yes 


output 


















Integration: baseline too high 


Less abundant ions underrepre- 
sented 








No 


UN 


UN 


UN 


Integration: baseline too low 


Less abundant ions overrepresented 






+ + + 


No 


UN 


UN 


UN 


"Abundance sensitivity" 


Underrepresentation of more abun- 
dant ions 








Yes 


UN 


Yes 


UN 



TOF, time-of-flight mass analyzer; sector, magnetic sector instrument; NA, not applicable; +, increases (slightly); + , increases 

(moderately); + ++, increases (strongly); -, decreases (slightly); --, decreases (moderately); , decreases (strongly); Var, variable; UN, 

unknown. 



nonlinearity might be correctable by use of suitable 
standards, if these could be synthesized, or by use of 
improved multiplier-detectors. Detector nonlinearity 
would also predispose to abundance-sensitivity effects: 
if abundances of different isotopomers span a large 
dynamic range, then the slope of detector output vs. 
injected material will be different for each. A third 
possibility is the occurrence of chemistry in the ion 
source (10, 11, 30). Chemical reactions such as hydro- 
gen abstraction or addition could alter mass isotopomer 
abundances. This problem could be addressed by using 
derivatives that are unlikely to undergo these reactions 
(such as fluorinated molecules lacking exchangeable 
hydrogens) or by use of different ionization conditions 
(e.g., metastable atom bombardment). This problem 
might also explain concentration effects on isotope 
ratios (increased concentrations in the ion source lead 
to more chemical interactions, as postulated in Refs. 
10. 11, and 30). 

Regardless of the physical mechanisms involved, 
some empirical observations may be helpful operation- 
ally. First, concentration effects (relative abundances of 
isotopomers varying as a function of the total amount of 
material injected onto the mass spectrometer) are a 
problem for most molecules and tend to be more 
pronounced the greater the dynamic range among the 
masses monitored. With glucose-pentaacetate, for ex- 
ample, in natural abundance samples the ratio of the 
Mo (0.8396) isotopomer to the Mj (0.1348) and Mz 
(0.0256) isotopomer abundances spans almost two or- 
ders of magnitude. In contrast, at higher values of 
precursor enrichment (e.g., p = 0.24), the ion envelope 
is more evenly spread out (Mq = 0.3060, M, = 0.4400 
and Mz = 0.2550, for a span of less than twofold), and 
observed concentration sensitivity is in fact less of a 
problem (Neese, R. A., R. Bandsma, and M. K. Heller- 
stein, unpublished observations). Second, the highest 
abundance isotopomers tend to be relatively underesti- 



mated as total ion abundance increases (as has also 
been observed with inorganic elemental analyses using 
isotope ratio mass spectrometers, Ref. 24). Third, and 
counterintuitively, masses with very low baseline iso- 
tope abundances (i.e., higher masses with few or no 
natural isotope abundances) can be analj^ically unde- 
sirable when they have to be compared quantitatively 
with higher abundance masses. Although it may seem 
attractive to avoid having to subtract baseline values, 
the extremely large dynamic range maximizes the 
relative concentration sensitivity of different isotopom- 
ers and may thereby lead to concentration-sensitivity 
for isotopomer abundances. 

In summary, a number of analytic factors influence 
the relative abundances of mass isotopomers measured 
in an envelope (Table 2). These factors range from 
chemical events in the ion source (abstraction of H) to 
performance of the mass analyzer (ion transmission 
efficiency, mass resolution, peak tailing), characteris- 
tics of the ion detector (velocity dependence of multi- 
plier, nonlinearity due to threshold or saturation ef- 
fects) , or accuracy of the integration software (baseline 
value that is subtracted, peak-fitting algorithm used). 
Until basic mass spectrometry research identifies the 
causes of quantitative inaccuracy, it will not be possible 
to correct post hoc in a definitive manner for substan- 
tial deviations from expected isotopomer abundances. 
Other strategies are therefore required. The most 
important of these in practice involve assessment of 
instrument performance and analytic techniques for 
prevention of inaccuracy. 

Strategies for Evaluating Quantitative Instrument 
Performance and Data Acceptability 

In the absence of methods for salvaging inaccurate 
analyses, the practitioner can establish criteria for 
acceptable accuracy and reject data that fail to meet 
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these criteria. Two general approaches can be used for 
evaluating instrument accuracy and data acceptability. 

Measurements on natural abundance standards. The 
theoretically expected mass isotopomer abundances 
can be calculated for any natural abundance molecule 
or ion fragment of known chemical composition (see 
appendix) . The impact of variations in natural isotope 
abundances is extremely small (Fig. 2) , so that natural 
abundance distributions should reliably predict mea- 
surements from standards. The simplest test of instru- 
ment accuracy is therefore how closely the measured 
isotopomer abundances in a natural abundance sample 
match the theoretically expected values. Minimum 
accuracy criteria for data acceptability can be empiri- 
cally established on the basis of statistical consider- 
ations (e.g., all isotopomers must be within 2% of their 
theoretical natural abundance values for the data to be 
included; thus 0.0025 for an My of palmitate-methyl 
ester or 0.0005 from Mz of palmitate-methyl-ester or 
glucose pentaacetate), taking into account the degree of 
accuracy required for a particular application, which in 
turn is influenced by the isotopomer enrichments 
achieved in the experiment. Obviously, accuracy needs 
to be at its best if the actual isotope enrichments 
present are low. One can formalize these criteria by 
simulating the effect of different types of errors, e.g., 
constant fraction or constant amount errors on calcu- 
lated parameters (see Achievement of adequate enrich- 
ments of masses ofinteresi) . 

Comparison of baseline measurements to theoretical 
abundances represents a simple strategy for establish- 
ing whether mass resolution, nonlincarity of detector 
response, or ion chemistry is a significant problem in a 
particular analysis. The main caveat is that analysis of 
natural abundance molecules cannot ensure that linear- 
ity is maintained as the relative abundances of differ- 
ent ions change, i.e., that accuracy will be maintained 
for labeled, biological samples. The latter may require 
use of labeled standards. Also, it should be noted that 
performance for mass isotopomers with essentially zero 
natural abundance values (e.g., M4 glucose) cannot be 
evaluated by this method, which makes it more difficult 
to assess quality of data if these masses are subse- 
quently monitored in labeled molecules (32). 

Use of isotopically enriched mass standards to estab- 
lish linearity of response for different isotopomers. One 
can use labeled standards to assess linearity and 
quantitative accuracy at various masses. [l-^HJ- and 
[6,6-2H2]glucose or [l-'^C]- and [l,2-i3C2]palmitate, for 
example, can be purchased and mixed with natural 
abundance molecules to generate mass standards. Af- 
ter correction for the isotopomeric envelopes actually 
present in labeled molecules (see appendix), standard 
curves representing different values of p and /can be 
analyzed. Problems with this approach include the fact 
that standards are themselves not 100% isotopically 
enriched, so one has to use an instrumental measure- 
ment of enrichment at some point to evaluate instru- 
ment performance, which can lead to circularities of 
logic, and the fact that accuracy of pipetting, complete- 
ness and comparability of derivitization, and reproduc- 



ibility of the injector can influence the standard curves 
observed. Moreover, this approach does not provide a 
definitive method of correcting for inaccuracies of abun- 
dance measurements, but only of evaluating whether 
inaccuracy is present and severe enough to invalidate 
results. Nevertheless, this approach is useful for evalu- 
ating quantitative linearity over a range of isotopomer 
abundances (3, 28) and is particularly useful when a 
mass isotopomer being monitored has insufficient natu- 
ral abundances to be measurable in baseline samples. 

Guidelines for Optimizing Quantitative Mass 
Spectrometric Analyses 

There are some useful experimental guidelines that 
have proven helpful in preventing instrument inaccu- 
racy and reducing the impact of instrument perfor- 
mance on kinetic parameters with the use of MIDA 
(26, 30). 

Attention to concentration effects. Publications that 
use MIDA should show the concordance of baseline 
isotopomer ratios with expected abundances (or show 
some other index of accuracy) and should demonstrate 
that concentration effects were considered and avoided. 
Failure to report these data makes technical evaluation 
of results difficult or impossible. 

Achievement of adequate enrichments of masses of 
interest. A second factor that is to some extent under 
the control of the investigator is the enrichments 
achieved in mass isotopomers of interest. Improving 
enrichment relative to background abundances im- 
proves reliability of parameter estimates and reduces 
sensitivity to analytic error, as we have discussed 
previously (26, 37). As an empirical rule, any isotopic 
enrichment <0.0050 (0.50 MPE) is problematic for 
MIDA calculations because of the large coefficient of 
variation that is unavoidable at such low enrichments 
(e.g., imprecision or inaccuracy of even ±0.0010- 
0.0020 represents a >20-40% coefficient of variation) . 
A formal error analysis at 0.0050 M2 enrichment in a 
sample polymer (palmitate) is shown (Fig. 7). These 
guidelines are dependent on the state of the technology, 
of course; if mass spectrometers improved accuracy and 
precision by an order of magnitude, the lower limit of 
acceptable incorporation would change accordingly. 

The investigator can both estimate and influence the 
enrichments that are attained in an experiment (Fig. 
8). The determinants of isotopomer enrichments in a 
polymer are n (the number of subunits), p (precursor 
pool labeling), and /"(fractional synthesis). The value 
for /is generally one of the unknowns being investi- 
gated, so this cannot be manipulated by the investiga- 
tor. The value for p can be controlled, however, by the 
dose of isotopic substrate administered. Here the bal- 
ance is between avoiding an undesired substrate (non- 
tracer) effect in the pool of interest vs. achieving 
optimal enrichments in the product. It is sometimes 
possible to select the value of n that will be analyzed 
(e.g., for peptide fragments of a protein into which 
leucine or another labeled amino acid has been incorpo- 
rated, Ref. 28), but this is more often not adjustable. 
The general relationship between p, n, and enrich- 
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Fig. 7. Simulated error-analysis at low enrichment of an isotopomer 
In a polymer, in this case, Mz of methyl-palmltate at 0.0050 isotopic 
enrichment. Effects on calculated parameters (p and of an anal3ftlc 
error of ±0.0010 in LAz are simulated across the range of true pfrom 
0.00 to 0.10 (at steps of 0.005). Calculated values of p and /(calc. p, 
calc. fj are shown as a percentage of correct values of p and / 
Different models of analytic error are compared: e.g., misidentifica- 
tion of Mz vs. Mi isotopomers {•-•); altered linearity of detector 
resulting in quantitative error in Mz abundances only (O-O). Solid 
lines, added to Mz; dashed lines, subtracted from Mz. 




Fig. 8. Simulated relationship between n (polymer length as number 
of subunits), p, and relative abundance of double-labeled species in 
the case of a general polymer (for simplicity, natural abundance 
background is not included in the calculations). 



ments of double-labeled species is shown in Fig. 8. 
Simulations (see appendix) can provide useful informa- 
tion about the optimal dose of isotope to administer. 

Use of standard curves. Most non-MIDA kinetic 
applications in stable isotope-mass spectrometry (e.g., 
dilution measurements of an infused tracer to deter- 
mine metabolite flux) do not ultimately require instru- 
ment accuracy. Standard curves of the analyte can be 
made simply by mixing known amounts of labeled and 
unlabeled molecules in different proportions; the stan- 
dard curve allows measured isotope enrichment to be 
converted to true proportions of labeled species present 
(tracer/tracee). The problem for biosynthetic measure- 
ments is that a simple standard curve cannot be easily 
made from available reagents, because there is no 
"standard" molecule for comparison. Each biological 
experiment results in a unique combination of p and f, 
thus a unique pattern of mass isotopomers are present 
in the population of molecules. At a minimum, one 
would have to establish the performance of the instru- 
ment by creating a standard curve of p vs. measured ion 
abundances. It would be better to mix enriched mol- 
ecules with unlabeled molecules to simulate end- 
product dilution (/). Testing of three-dimensional stan- 
dard curves in this manner requires substantial effort 
and analytic time, however, and does not correct for 
concentration effects (without imposition of another 
analytic dimension for the standards). Although mix- 
tures of labeled standards can be useful for assessing 
instrument performance (3, 28 and previous discus- 
sion), the role of higher-dimensional standard curves in 
MIDA applications remains uncertain. 

Does the MIDA Calculation Algorithm Used Affect 
the Final Parameter Estimates, in Practice? 

As discussed above, the published algorithms for 
calculating biosynthetic parameters by MIDA from 
experimental data are similar in many respects and 
give essentially identical parameter estimates in the 
theoretical case (5). One potential difference is the 
capacity of different calculation approaches to adjust to 
instrument inaccuracy. Some calculation algorithms 
subtract measured baseline (natural abundance) data 
from enriched sample data before transforming sample 
data or fitting them to a model; other algorithms do not 
subtract measured baselines but correct for theoretical 
baselines or simply use the enriched sample data as 
they stand. In the former case, the goal is to account for 
instrument performance on a sample-to-sample basis; 
in the latter case, ideal instrument performance is 
assumed. The capacity of these two strategies to adapt 
to instrument inaccuracy can be simulated (Table 3). In 
this simulation, we considered several potential types 
of measurement inaccuracy: spillover between adjacent 
ions, constant contaminant at a particular mass-to- 
charge ratio {m/z), nonlinearity of the detector at 
different m/z values. The impact on parameter esti- 
mates of using theoretical values vs. measured base- 
lines was modeled. It is clear that use of measured 
natural abundance values in all cases reduces the error 
in parameter estimates, although not completely. This 
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Table 3 . Simulated analytic errors in glucose 
pentaacetate: effects of baseline values used 
on calculated values off 



Type of Error Involved 


BSLN 
Used 


Calculated / 

% actual value) 


5% Spillover (Ml to Mo) 


Theoretical 


80 




"Measured" 


88 


5%Spmover(MotoMi) 


Theoretical 


126 


"Measured" 


114 


Constant contamination (0.0050 in Mi) Theoretical 


107 




"Measured" 


100 



Simulations were carried out by Imposing a unidirectional constant 
fraction error (5% of Mi into the Mq isotopomer or 5% of Mq into the 
Ml isotopomer) or a constant value error (0.0050 of the total envelope 
into the Mi isotopomer) on labeled glucose pentaacetate (p = 0.15). 
Ai4i, Ai4z, and parameter estimates pand /were then determined, 
after subtraction of either theoretical baseline (BSLN) values for 
glucose pentaacetate ions monitored [mass-to-charge ratio (miz) 331 
to 333, see text] or "measured" BSLN values (i.e., imposing the same 
Instrumental error as for labeled glucose). 

conclusion was also reported by Chinkes et al. (5). Thus 
baseline corrections cannot substitute for accurate 
quantitation by mass spectrometers but may reduce 
error in the final calculated values. 

Is There Experimental Evidence Supporting 
or Contradicting the MIDA Model? 

A general measurement technique such as MIDA 
needs to be tested under controlled conditions both in 
vitro and in vivo. Strong experimental evidence in 
support of the method has been generated in vitro. Lee 
et al. (22) synthesized glucose-pentaacetate from 
['^C] acetic-anhydride at known enrichments and then 
imposed dilutions by known amounts of unlabeled 
glucose-pentaacetate. The calculated values of p and f 
were nearly identical to the expected values. We (28) 
performed an analogous test by synthesizing an oligo- 
peptide (SVVLLLR) from [^Hs] leucine at known enrich- 
ments and imposing subsequent dilutions by unlabeled 
peptide. Again, observed values were extremely close to 
predicted values. The combinatorial probability-mass 
isotopomer quantification model is clearly an accurate 
theoretical description of polymerization biosynthesis 
under controlled conditions. 

Testing a method definitively in vivo is not as straight- 
forward. Comparisons with alternate methods, e.g., the 
close similarity of lipogenic estimates by 2H2O and 
MIDA (reviewed in Ref. 16) or of cholesterogenesis by 
sterol balance and MIDA (25) are not definitive, be- 
cause the comparison methods may be flawed as well. A 
useful result was presented by Kelleher et al. (20): 
cholesterogenesis approached the 100% value expected 
in exponentially growing cells in culture. 

An important issue to discuss is the proper interpre- 
tation of apparently dissonant results with an in vivo 
method like MIDA. A recent example is the intriguing 
observation of Aarsland et al. (1) that total hepatic 
lipogenesis (based on MIDA on circulating very low 
density-triglyceride fatty acids) was less than l/20th of 
the rate of net whole body lipogenesis (based on indirect 
calorimetry) during massive carbohydrate overfeeding 



in humans. Because net lipogenesis (synthesis minus 
oxidation) cannot be greater than unidirectional lipogen- 
esis, is this evidence that MIDA failed to give a 
physiologically possible answer? This conclusion would 
be incorrect: the authors (1) instead concluded that 
lipogenesis occurred in tissues or sites not immediately 
communicating with circulating very low density lipo- 
protein triglycerides during an 8-h isotope infusion. 
Thus lipogenesis in adipose tissue or lipogenesis enter- 
ing the hepatic cytosolic storage pool could account for 
the unmeasured lipogenesis in the whole body. Physi- 
ological explanations must be excluded when the valid- 
ity of MIDA is evaluated in vivo. 

Analytic factors must also be excluded, particularly 
those related to instrument accuracy and concentration 
effects on measured isotope abundances (10, 11, 30). In 
our view (7, 26), analytic and experimental design 
factors may to some extent explain recent discordant 
results (32) of the use of MIDA for gluconeogenesis. 
Although this particular question has not been re- 
solved, in general it is essential for manuscripts testing 
MIDA to demonstrate accuracy, avoid abundance- 
sensitivity effects, and explicitly describe the measures 
used to do so. Claims that a biosynthetic system is not 
adequately described by a combinatorial model need to 
be carefully evaluated using criteria other than just 
dissonant final estimates. 

SOME FUTURE DIRECTIONS FOR MIDA 

A number of future directions can be considered for 
combinatorial probability techniques. Applications with 
large, heteronuclear polymers (e.g., proteins, polynucleo- 
tides) represent a challenge but potentially include 
some of the most important molecules in biology. Spe- 
cial problems are introduced for heteronuclear poly- 
mers, such as proteins, for example, because their large 
size and composition (including 20 or more amino acid 
subunits) require the investigator to identify repeats of 
a particular homonuclear subunit and then isolate this 
subunit from the intact molecule (28). Although tech- 
niques such as electrospray ionization are able to 
introduce intact proteins into the gas phase for mass 
spectrometry, resolution of individual mass isotopom- 
ers is not generally possible with most mass analyzers, 
especially for the multiply charged species generated 
by electrospray ionization. We (28) have successfully 
used the strategy of enzymatically hydrolyzing proteins 
to proteolytic fragments, and then collecting a selected 
fragment that includes a homonuclear stretch (e.g., 3-5 
leucines out of a 7- to 15-amino acid stretch). An 
analogous approach might be applicable to DNA or 
RNA samples, which have the advantage of containing 
only four different subunits. 

Another potentially powerful direction for MIDA may 
be to exploit information about p, not just f Establish- 
ing the timing of biosynthetic events (e.g., during 
embryonic development) might be possible by finger- 
printing polymers according to their precursor pool 
enrichment after generating a time gradient within the 
precursor pool. Measuring the tissue location of a 
biosynthetic event (e.g., whether cholesterol in high- 
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density lipoprotein was synthesized in peripheral tis- 
sues or the liver, representing reverse or forward 
cholesterol transport) might also be possible by analyz- 
ing nonoverlapping values of p within an isotopomeric 
envelope, if label can be delivered at different rates to 
different tissues. Testing or correcting for labeling 
gradients within a tissue may also be possible by 
analyzing different portions of an isotopomeric enve- 
lope for divergent values of p (Fig. 5Q. 

Conclusion 

An algorithm (see appendix) for calculating fractional 
abundances of mass isotopomers in complex mixtures 
of labeled and natural abundance polymers can be used 
to construct reference tables that allow inference of the 
biosynthetic parameters p and f. The consequences of 
deviations from various assumptions of the MIDA 
model and the practical limits of the technique can also 
be tested in this manner. MIDA appears to be robust in 
the face of a number of deviations from its central 
assumptions. The combinatorial-probability approach 
imposes minor or no constraints with regard to the 
isotope-labeled substrate that can be administered, the 
length of the polymer that can be studied, the presence 
of contaminating M-i fragment ions, the existence of 
inhomogeneity or isotopic disequilibrium in the mono- 
meric precursor pool, the possibility of variations in 
natural ^^C abundances, or the existence of isotope 
discrimination. More problematic is analytic inaccu- 
racy, which can be partially diagnosed by incorrect 
baseline isotopomer abundances and which is often 
exacerbated by concentration-sensitivity. Whether ana- 
lytic imprecision or deviations from the assumptions of 
the simple MIDA model have an impact on biosynthetic 
parameters in a particular case can be tested by 
individualized simulation. 

APPENDIX: GENERAL CALCULATION ALGORITHM 
FOR MEASURING POLYMERIZATION BIOSYNTHESIS 
BY MIDA2 

Our previous descriptions of the MIDA method (13, 14, 16, 
17) provided the underlying theoretical basis but presented 
only an overview of the mathematical approach. The actual 
computer-based calculation algorithm that we have used for 
the past several years has not been described since our 
original presentation of the technique. We will describe a 
systematic and general calculation algorithm that can be 
used with a personal computer to generate expected frac- 
tional abundances of mass isotopomers for polymers contain- 
ing labeled subunits; we will then show how these theoretical 
fractional abundances can be made into a "reference table" 
from which one can convert experimental data on mixtures of 
natural abundance and isotopically enriched polymers into 
biosynthetic parameters (pand /); and finally, we will demon- 
strate how to perform a MIDA biosynthetic calculation from 
experimental data. How the algorithm can be used to evalu- 
ate limitations of the technique or deviations from the model. 



or to correct for certain deviations from the model, is de- 
scribed in the text of this article. 

Calculation Algorithm 

We first present a formal mathematical algorithm for 
calculating the fractional abundances of mass isotopomers 
resulting from mixing natural abundance molecules with 
molecules newly synthesized from a pool of labeled monomers 
characterized by the parameter p. A mixture of this type can 
be fully characterized by f, the fraction new, and p. This 
algorithm is presented in a step-wise fashion, beginning with 
the simplest calculation, a molecule synthesized from a single 
element containing isotopes with the same fractional abun- 
dances that occur In nature and not mixed with any other 
molecules. We then proceed to molecules containing more 
than one element with all isotopes at natural abundance; 
then to nonpolymeric molecules containing different ele- 
ments, some of which are in groups whose isotope composi- 
tion is not restricted to natural abundance but is variable; 
then to polymeric molecules containing combinations of repeat- 
ing chemical units (monomers), wherein the monomers are 
either unlabeled (containing a natural abundance distribu- 
tion of isotopes) or potentially labeled (containing an isotopi- 
cally perturbed element group); and finally to mixtures of 
polymeric molecules, composed of both natural abundance 
polymers and potentially labeled polymers, the latter contain- 
ing combinations of natural abundance and isotopically per- 
turbed units. The last-named calculation addresses the condi- 
tion generally present in a biological system, wherein polymers 
newly synthesized during the period of an isotope incorpora- 
tion experiment are present along with preexisting natural 
abundance polymers, and the investigator is interested in 
determining the proportion of each that is present to infer 
synthesis rates or related parameters. 

Molecules containing only a single element at the fractional 
isotope abundances that occur in nature. A first step is to 
calculate the isotope pattern of a hypothetical molecule 
composed of the element / with the fractional abundances 
0/= 0.989, '/= 0.01, 2/= 0.001. In the limit the molecule is a 
single atom, and the isotope pattern in a collection of element 
/ is obtained directly from the fractional abundances of the 
isotopes: Aq = 0.989, Ai = 0.01, and A^ = 0.001. However, for a 
molecule composed of N I atoms (/n), the mass isotopomer 
distribution (isotope pattern) is obtained from the multino- 
mial distribution. This distribution is a function of two 
variables; N (the number of / atoms) and the fractional 
abundances of the isotopes, which in this case are not 
perturbed from natural abundances. The contribution at each 
mass unit is obtained by individually summing the fractional 
abundances of all the isotopologues with that particular 
nominal mass. In this fashion one obtains the fractional 
abundance that would be observed at that mass with a mass 
spectrometer operating at unit mass resolution, arising from 
the conglomerate of isotopologues with that nominal mass. 
The frequency (Fj of any given isotopologue in the summed 
total abundance of all of the isotopologues is given by the 
multinomial distribution [as adapted from CRC Standard 
Mathematical Tables (26th ed., pg 519); see Eg. Al]. For 
example, the probability of observing the Ig molecule com- 
posed of °/2'/5^/i is calculated as follows 



2 Richard A. Neese, Dennis Faix, Kennetii Caldwell, and Marc K. 
Hellerstein were the authors of the appendix. Dr. Caldwell is 
deceased as of August 1997. 
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Note that this is the fractional abundance for one particular 
isotopologue of nominal mass m^ due to the isotopologue 
"Iz^Is^I], not the fractional abundance of the particular mass 
isotopomer (M?) for the molecule, since there are other 
isotopologues that have a nominal mass of 7777 (e.g., ^I^^I^I^. 
Thus the multinomial distribution allows the calculation of 
the probability or fractional abundance of any unique combi- 
nation of isotopes (any isotopologue) in a molecule composed 
of a single element. To obtain the final probability of any 
particular mass isotopomer, one must sum the probabilities 
given by the multinomial distribution for each of the possible 
isotope combinations of the element. We note here that /s 
would become an element group (defined as a group of atoms 
in a molecule attributed to a given element and sharing the 
same fractional abundances of the isotopes of the element; the 
fractional abundances of isotopes of the element, whether 
natural or perturbed, are by definition the same for all atoms 
in an element group, and a molecule may contain more than 
one element group for a given element) . This element group is 
part of a molecule if that larger molecule contained eight / 
atoms and each of the / atoms were from the same isotopic 
pool. 

Molecules containing more than one element at natural 
isotope abundances. Next, suppose that probability distribu- 
tion calculations have been done for the element groups C2, 
H4, and O2 as if they were isolated molecules, each composed 
of a single element, as described in the previous section. How 
then is the mass isotopomer distribution for C2H4O2 (acetic 
acid) obtained? One way is to multiply the fractional abun- 
dances of the components of each element group together and 
then to sum the fractional abundances of the specific isotopo- 
logues according to their nominal mass. Accordingly, to obtain 
the Ml components of acetic acid, the M component of each 
element group is combined with the Mq components of the 
others: (i^C'^C, IH4. i^Oz), ('^Cz, 'Hs^H, ^^Oz), and (12C2, 'H4, 
'^O'^O). Each of these unique combinations has a probability 
of occurring that can be calculated by using the multinomial 
distribution {Eq. Al); they are multiplied together to obtain 
the fractional abundances of the isotopologues (see Eq. A2). 
The fractional abundance (A{) of the Mi mass isotopomer of 
acetic acid is then the sum of these fractional abundances, 
representing the sum of all the isotopologue fractional abun- 
dances {Eq. A2d). 

Although calculating the fractional abundances of indi- 
vidual isotopologues and summing them according to mass is 
a reasonable approach for obtaining a mass isotopomer 
distribution in simple molecules, a systematic, step-by-step 
algorithm is needed for more complex molecules. We now 
describe an algorithm, first in general terms and then by 
using acetic acid as an example. The strategy, alluded to in 
the previous example, is to disaggregate the molecule concep- 
tually into element groups (e.g., C's, H's, O's), calculate the 
mass isotopomer distributions for each element group, and 
then multiply the appropriate components of each element 
group together in a systematic fashion to attain the final 



mass isotopomer distribution of the molecule. A computer 
algorithm can do this by exhaustively combining all of the 
components of the first two element groups until the mass 
isotopomer distribution corresponding to the hypothetical 
molecule composed of those two element groups is obtained, 
and then by recording these abundances in an array. The 
process is then repeated for the next element group by 
combining its fractional abundances with those in the array. 
The resulting fractional abundances in the new array are 
then combined with the fractional abundances of the compo- 
nents of the next element group, and so on. A key aspect of 
this procedure is specifying all the combinations that contrib- 
ute to the fractional abundances at each nominal mass during 
each iteration of the algorithm. For example, only Mq values 
from the array and the element group contribute to the Mq 
values in the new array, so they are multiplied together; 
however. Mi and values contribute to the new Mi, so they 
must be multiplied each time. Likewise, Mq, Mi, and M2 
values from the two element groups contribute to the new Mz 
values, so Mq values are cross multiplied with M2 values, and 
Ml values are multiplied together. When all of the element 
groups have been combined together, the final distribution of 
fractional abundances for the whole molecule is attained. In 
the computer subroutine, this multiplication and summing 
process is carried out by use of loops, one nested inside the 
other. Because each addition of an element group increases 
the number of mass isotopomers and because the fractional 
abundances of mass isotopomers greater than approximately 
Mio are typically extremely low, an arbitrary limit may be set 
so that abundance calculations are not performed beyond this 
or some other set limit. 

The example of acetic acid demonstrates how the computer 
algorithm is used (Table 4). The two element group distribu- 
tions C2 and H4 are combined to obtain the mass isotopomer 
distribution of the intermediate C2H4 molecule. This new 
distribution is in turn combined with the O2 distribution to 
form the final mass isotopomer distribution for acetic acid. 
The Mo value for C2H4 (0.97771) is multiplied by the Mq value 
for O2 (0.99519) to yield the fractional abundance of the Mo 
mass isotopomer for acetic acid [Aq = 0.97300). Because only 
Mo mass isotopomers of the element groups contribute to that 
mass isotopomer for the final molecule, the calculations are 
complete for A^, and the routine moves on to the next one, the 
Ml mass isotopomer. In this manner A^ to A^ are calculated 
(Table 4). 

To facilitate combining two distributions, a single mass 
isotopomer distribution can be thought of as a set of values 
corresponding to the fractional abundances of each of the 
possible mass isotopomers of the molecule; any distribution 
(for an element group or molecule) is a set with n terms. 
Mo— Mn--i. Two distributions, A and B, will have n^ and ns 
terms, respectively, whereas a distribution resulting from the 
combination of A and B will have Hab terms, where n^a = 
Ha+hb - 1. The formula for combining two distributions, A 
and B, into a new single distribution AB is given in Eq. A3. 



Table 4 . Mass isotopomer abundance distribution calculations for natural abundance acetic acid ( C2H4O2) 

Mass Isotopomer (M,) 



tsiemeni: 


Mo 


Ml 


Mj 


M3 


Mt Mi Me 


Mj 




C2 


0.97832 


0.02156 


0.00012 










H4 


0.99938 


0.00062 


0.00000015 


1.5x10-11 


5.9x10-16 






O2 


0.995186 


0.000738 


0.00407 


0.00000151 


0.000000416 






C2H4 


0.97771 


0.02216 


0.00013 


7.7X10-« 


1.8X10-11 1.8x10-1= 






C2H4O2 


0.97300 


0.02277 


0.00413 


0.00009 


4.6x10-6 9.3x10-8 5.5x10-10 


3.2 X 10-" 


7.4x10-1' 
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Nonpolymeric molecules containing an element group in 
which isotope composition is perturbed from natural abun- 
dances and is variable. Next, the calculation of mass isoto- 
pomer distributions is shown for nonpolymeric molecules that 
contain an element group enriched in a particular isotope, 
that is, perturbed from the natural fractional abundances of 
its isotopes. For a molecule I^y, where a number {w) of the / 
element group is enriched in the calculation is nearly the 
same as described already except that the w atoms of /must 
be considered as a separate element group (*I) from the other 
(x- w) atoms of /in the molecule. The molecule then becomes 
*Iy,Ix-wJy with the appropriate distribution calculated for 
each element group and the element group distributions 
combined to generate the final distribution of fractional 
abundances of mass isotopomers in the whole molecule. The 
first step operationally, as described in the previous ex- 
amples, is to break the molecule into element groups, calculat- 
ing Jy and 4- w from the natural isotope distributions of 
elements J and /. It is important to emphasize that Jy and 
4- will have constant mass isotopomer abundance distribu- 
tions even as the enriched elemental group (*/) changes. This 
invariant portion will be referred to here as the constant or 
invariant mass isotopomer abundance distribution of the 
molecule, and it can be treated as though it were a chemical 
derivatizing agent attached to the variable moiety. Dividing 
the isotopically perturbed molecule into constant and vari- 
able moities greatly simplifies the calculation algorithm and 
is the motivation for dividing the elements in a molecule into 
element groups (one of which may be constant while the other 
is variable). 

After the constant distribution is calculated, the distribu- 
tion in the variable moiety is calculated. It is useful to divide 
the */w element group itself into two populations. The first 
population is the proportion (p) of the element group that is 
the perturbed portion, made up of isotopically labeled atoms 
(e.g., '/). The second population is the proportion 1 - p of the 
element group that is made up of the natural abundance 
atoms. The abundance distribution of the labeled proportion 
('/„) is simple^: ^4^ = 1.0. For the natural abundance propor- 
tion of the element group ("^/w), the distribution is calculated 
in the manner described above, by using the multinomial 
distribution (Eq. Al) with natural abundance values for w I 
atoms. Because these two populations are physically mixed at 
the site of synthesis of the molecule from the combined 
distribution must be calculated. This combined distribution is 
calculated by correcting for the proportion of each present at a 
given value of p: the labeled distribution is multiplied by p 
and the natural abundance distribution by 1 - p, and the 
fractional abundances of the corresponding nominal masses 
are sunrimed. Thus a weighted distribution is generated for 
*/w. This new distribution is the fractional abundance distri- 
bution for at that particular value of p and varies as a 
function of p. The */„ distributions over a broad range of 
values for p are then combined with the constant distribution 
to calculate the final molecular distribution as a function of p. 

By way of illustration, the mass isotopomer fractional 
abundance distribution is calculated for a theoretical pool of 
[^Ha] leucine molecules, which has been synthesized chemi- 
cally. The element groups for these leucine molecules 
(C6H13O2N) are as follows: Ce, Hjo, O2, and N (the constant 
portion, all at natural isotope abundances) and *H3 (the 
perturbed or variable portion, with a variable distribution). 



3 We will assume 100% isotopic labeling by +i/in the material. If 
the material is <100% the actual distribution of i/„,is calculated 
by using Eq. A3. 



The Ce, Hio, O2, and N distributions are calculated individu- 
ally and combined [Eq. A3) to form the constant distribution, 
which then must be combined with the distribution of the 
variable portion of the molecule. For the variable portion, the 
fractional abundance of ^Hs in *H3 varied from zero (i.e., all 
*H3 values are at natural abundance) up to a reasonable point 
of interest, when the fractional abundances of Mq— M3 are 
calculated for each value of p. For example, if p = 0.04 (4% of 
the *H3 element group consists of ^Ha), 96% will have the 
natural abundance distribution for *H3 {Aq = 0.99953, Ax = 
0.00047, Ai = 0.000, A^ = 0.000), and 4% will have a 
distribution for *H3 of A3 = 1.0. Combining these two 
distributions by using a weighted average, or linear combina- 
tions, gives a final *H3 distribution at p = 0.04 of Aq = 
0.95955, Al = 0.00045, Ai = 0.00, and A3 = 0.04. This 
distribution is then combined with the constant distribution 
to obtain the final molecular isotopomer distribution, as seen 
in Table 5 for p = 0.04. Each *H3 distribution at a given value 
of p is combined with the constant distribution to give a 
distribution of mass isotopomer abundances for the molecule 
leucine as a function of p CTable 5). 

Polymeric molecules containing repeating monomeric sub- 
units, with some of the monomeric subunits composed exclu- 
sively of element groups at natural isotope abundances while 
other monomeric subunits contain an isotopically perturbed 
element group. The next step is to calculate the fractional 
abundances of mass isotopomers in polymers containing 
combinations of monomers, which themselves are either 
unlabeled (containing only element groups at natural isotope 
abundances) or labeled (containing an isotopically perturbed 
element group). The calculation procedure is an extension of 
the principles described so far. A general discussion is pre- 
sented first, followed by an example of a peptide containing 
repeated subunits of leucine, .either at natural abundance or 
containing PH3] leucine (see above). 

Consider a polymer composed of /^subunits, i.e., (IxJ^r If 
a number [w) of / atoms (w<x) in each I^y subunit are 
enriched in +'/(*/), whereas all other atoms are at natural 
isotope abundances (/and J], the calculation for each subunit 
I^Jy is as just described for a nonpolymeric molecule that 
contains an element group enriched in a particular isotope. 
The mass isotopomers of the polymer {IxJy)z then can be 
calculated and will depend upon p. Several mathematically 
interchangeable algorithms exist for calculating the distribu- 
tion of mass isotopomer abundances for a polymer of this 
type. The calculation approach that we believe to be the 
simplest computationally is to treat the isotopically per- 
turbed element group in the polymer as a single, discrete 
(inseparable) unit, with its own mass isotopomer abundance 
distribution that is then combined with the constant distribu- 
tion of the remainder of the polymer. The first step, as with 
the described nonpolymeric leucine example, is to break the 
polymer into element groups. The / positions that can be 
enriched are treated as a separate polymer (*/Jz. This is 
convenient because the mass isotopomer abundance distribu- 
tion of (*/Jz will vary as a function of p, whereas the 
remaining element groups Jy^ and remain at natural 

isotope abundances and therefore have a constant abundance 
distribution of mass isotopomers, even as the percentage of 
isotopically perturbed subunits changes. This constant or 
invariant distribution of the poljmier is treated identically as 
in simple, nonpolymeric molecules (Table 4). Therefore, one 
can first calculate the mass isotopomer abundance distribu- 
tions for the element groups that are invariant. These are 
then combined, one at a time, until the constant distribution 
is calculated. 
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Table 5 . Mass isotopomer fractional abundance distributions for leucine pool at varying proportions 

(p) of labeled pHsJleucine 



Fractional Abundance 



Mass Isotopomer 


M) 


M, 




Mj 


M4 


Ms 


Me 


Nominal mass (/j^z) 


131 


132 


133 


134 


135 


136 


137 


Exact mass 


131.10 


132.10 


133.10 


134.10 


135.10 


136.10 


137.10 


P 

0.00 


0.92657 


0.06723 


0.00588 


0.00031 


0.00001 


0.00000 


0.00000 


0.02 


0.90804 


0.06588 


0.00577 


0.01884 


0.00135 


0.00012 


0.00000 


0.04 


0.88950 


0.06454 


0.00565 


0.03738 


0.00269 


0.00024 


0.00001 


0.06 


0.87097 


0.06319 


0.00553 


0.05591 


0.00402 


0.00035 


0.00002 


0.08 


0.85244 


0.06185 


0.00541 


0.07444 


0.00536 


0.00047 


0.00002 


0.10 


0.83391 


0.06050 


0.00530 


0.09298 


0.00669 


0.00059 


0.00003 



Calculations for the variable part of the polymer are 
more complicated. Its fractional abundance distribution will 
depend on the number of labeled subunits (*!„,) that happen 
to be incorporated. If we call the number of labeled subunits 
a, this element group may be composed entirely of labeled 
subunits (a = z), or it may have 1, 2, 3 and so on, labeled 

subunits (o- = 1, 2, 3 ), or it may have no labeled subunits 

(a = 0). The remaining unlabeled subunits (z-ct) contain only 
elements at natural isotope abundances. The general formula 
for this variable element group is therefore (*lJ)„(Qz-a- 

Accordingly, the next step is to calculate all the possible 
mass isotopomer distributions for a = 0 to z. If ct = 0, then the 
calculation is simply that of inserting the natural abundance 
values into the multinomial distribution to obtain a distribu- 
tion for wz / atoms (I,,^), as was calculated for the constant 
part of the polymer. If cr = 1, then the distribution is obtained 
by calculating the natural abundance distribution for vv(z-\) I 
atoms [Iw{^-\)] and the distribution for *I„ (Ai = 1.0) and 
combining these two distributions {Eq. E3) to generate the 
distribution of the variable portion of a polymer with a single 
labeled subunit. If ct = 2, the natural abundance distribution 
Uwi^z-z)] is combined with the distribution {Ai = 1.0) in the 
same manner. The process is repeated up to a = z. If all the 
subunits are labeled (ct = z), then the distribution for *I„^ 
becomes A^ = 1.0, and all other mass isotopomer fractional 
abundances are zero. Stated in general terms, when the 
isotopomer distribution for an element group containing a 
perturbed or labeled subunits is calculated, the distributions 
for (labeled) and Iw(z-^) (natural abundance) are calcu- 
lated and combined to give the distribution for that subunit 
combination in the element group. 

These distributions for each value of a serve as reference 
distributions and need only to be calculated once for any 
given polymeric molecule containing a particular labeled 
element group. They do not, however, by themselves reveal 
the actual distribution in the polymeric molecule, because the 
proportion of each of these theoretical species (ct = 0, 1 , 2, 3,...) 
actually present must also be calculated, and this is deter- 
mined by p, the fractional abundance of isotopically labeled 
subunits in the precursor pool. The multinomial, or, in this 
example, the binomial expansion {Eq. Al),\s used to calculate 
the distribution of subunits in the element group as a 
function of p and z, the total number of subunits in the 
element group {Eq. A4). Each value of p will be associated 
with unique proportions of a = 0, 1, 2, 3, ... subunits in the 
element group (Table 6A). 

Accordingly, one then multiplies the frequency of each 
subunit combination (i.e., the coefficient) by the mass isoto- 
pomer abundances characteristic of that combination (i.e., 
the distribution) to generate a final, weighted distribution. 
This is the (*/J^mass isotopomer abundance distribution for 



that value of p (Table QB). This distribution is combined with 
the distribution of the invariant moiety of the polymer (using 
Eq. A3) to give a molecular distribution for that value of p. 
The process is repeated for every pof interest {O^.xxxat steps 
.OQx) until a table of mass isotopomer fractional abundances 
of the polymer as a function of p is obtained (Table 7). During 
iterations as p is varied, only the probability of each subunit 
combination needs to be recalculated before multiplying, 
summing, and combining to form the fractional abundance 
distribution of mass isotopomers of the molecule at that value 
of p (Table 65). The reference table obtained (Table 7) reveals 
data for p vs. isotopomer abundances in the intact SWLLLR 
peptide. It should be emphasized that each molecule analyzed 
by MIDA requires its own reference table. This needs to be 
calculated only once, however, across the range of values for 
p; the table can then be used for all subsequent measure- 
ments on this molecule. Computational intensity is substan- 
tial but applies most often at the beginning of work on a 
biosynthetic problem. 

Mixtures of polymers. Finally, one can consider mixtures of 
polymers, as might be present during a physiological isotope 
incorporation experiment in vivo, wherein newly synthesized 
polymers are being added to a preexisting pool of polymers as 
part of the process of turnover. This calculation consists 
simply of generating linear combinations of the fractional 
abundances of each mass isotopomer.'' Thus one multiplies 
the fractional abundances of mass isotopomers in each poly- 
mer population times the proportion of the polymer pool 
represented by each population, and one sums the results. In 
this manner, the fractional abundances of each mass isoto- 
pomer in the mixture are weighted for the proportion of each 
polymer present. The net result is a linear combination of 
mass isotopomer abundances {Eq. AS). 

Generation of Reference Tables for Inferring Biosyntlietic 
Parameters from Experimental Data 

Rationale for use of ^fractional abundances in reference 
tables. As discussed in the text of this article and previously 
(13, 15, 20, 21), inference of the fractional biosynthetic 
contribution to a polymeric pool {!) requires knowledge of the 
proportion (p) of precursor subunits that were isotopically 
labeled after introduction of a labeled monomer. The math- 
ematical reason why p must be known to establish /is evident 
by graphical analysis (see accompanying manuscript. Fig. 
IQ. The universe of possible fractional abundances for any 
mass isotopomer in a polymer can be seen to depend on the 



This holds unless an incomplete ion spectrum is being monitored, 
in which case a correction must be applied to restore the molar 
combinations to linearity (see text and Eq. A9b). 
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Table 6. Mass isotopomer distributions in peptide SVVLLLR (singly charged, protonated) due to incorporation 
ofpH3]leucineatp = 0.10 







Combinations of natural 


abundance lei 








[^Hslleucine 
Units (a) 




omic Elemental 
Formula 


Monoisotopic 
Mass, Jim 


Fractional Contribution 
to Peptide (Coefficient) 




1 




C37N10O9H71 (2Ho) 
CsyNioOgHfiS PH3) 


799.6 
802.6 


(0.9)^(0.1)° = 0.729 
j^jx {0.9)2(0.1}' = 0.243 




2 




CstNioOsHbs PHe) 


805.6 


(^jx (0.9) '(0.1)^ = 0.027 




3 




C37N 


0O9H62 (^Hg) 


808.6 


j^)x(0.9)''(0. 1)3 = 0.001 




[^Hslleudne 


B. Fractional abundances of mass isotopomers M 
Mo (799) M(800) MjCSOI) A^s (802) 


to M,o (m/z 799 to 809) in peptide at p= 0.10 for PHJleucine 
M,(803) Af5(804) Me (805) M, (806) (807) M9(808) 


Mo (809) 


0 

Distribution 

Coefficient 

Subtotal 

1 


0.621963 

0.729 

0.453411 


0.285353 

0.729 

0.208022 


0.075466 

0.729 

0.055015 


0.014608 

0.729 

0.010649 


0.002273 0.000299 0.000034 
0.729 0.729 0.729 
0.001657 0.0000218 0.000025 


0.000003 0.000000 0.000000 
0.729 0.729 0.729 
0.000002 0.060000 0.000000 


0.000000 

0.729 

0.000000 


Distribution 

Coefficient 

Subtotal 

2 


0.000000 


0.000000 


0.000000 


0.622255 

0.243 

0.151208 


0.285195 0.075368 0.014579 
0.243 0.243 0.243 
0.069302 0.018314 0.003543 


0.002267 0.000298 0.000034 
0.243 0.243 0.243 
0.000551 0.000072 0.000008 


0.000003 

0.243 

8.4X10-'' 


Distribution 

Coefficient 

Subtotal 


0.000000 


0.000000 


0.000000 


0.000000 


0.622546 
0.027 

0.000000 0.000000 0.016809 


0.285037 0.075270 0.014551 
0.027 0.027 0.027 
0.007696 0.002032 0.000393 


0.002261 

0.027 

0.000061 


3 

Distribution 

Coefficient 

Subtotal 


0.000000 


0.000000 


0.000000 


0.000000 


0.000000 0,000000 0.000000 


0.622837 
0.001 

0.000000 0.000000 0.001024 


0.284879 
0.001 

0.000285 



Totals 

Fractional 0.453411 0.208022 0.055015 0.161857 0.070959 0.018532 0.020377 0.008249 0.002104 0.001024 0.000347 
abundance 

%Relative 100.00 45,88 12,13 35,70 15,65 4,09 4,49 1.82 0,46 0,23 0,08 
abundance 



Distribution of the invariant moiety of the molecule is corrected for the proportion of each combination of labeled to unlabeled leucine 
subunits (the coefficient), to generate abundances of each mass Isotopomer. 



two variables, p and f. By envisioning a plane drawn at the 
fractional abundance value for any (i.e., 0. 15, Fig. 1 C), it is 
apparent that an infinite number of solution pairs of p and / 
could result in any given value of A^. In a nonlinear system in 
which both p and /are unknown and variable, such as an in 
vivo biosynthetic system, it is therefore necessary to solve for 
p to linearize the relationship between fractional abundance 
Ax and / It is then possible to solve for f by simple algebraic 
means. 

The calculation problem that has to be solved for mixtures 
is that the polymers at natural isotope abundances also 
contribute to the abundances of mass isotopomers. There is, 
however, a simple way of correcting for natural abundance 
contributions and for linearizing the relationship between f 
and Ax (Eqs. A5-A7): when expressed as Afractional abun- 
dances, ratios among mass isotopomers are constant and 
independent of f {Eq. A7c). It follows that the isotopomer 
pattern of the population of newly synthesized molecules can 
be established by expressing the data as ratios among Afrac- 
tional abundances. The value for p can then be inferred 
directly from reference tables (see Table 7) , and the value for / 
is calculated easily by application of the precursor-product 
relationship. The calculation advantage of subtracting mea- 



sured baselines has also been recognized and incorporated 
into the algorithms of Lee et al. (22) and Chinkes et al. (5), 
although other published algorithms do not include a base- 
line correction (21). 

Another mathematical point concerning the use of Afrac- 
tional abundances is. worth mentioning. The sum of deple- 
tions (negative values) plus enrichments (positive values) 
compared with natural fractional abundances for all the mass 
Isotopomers analyzed must add up to zero, because fractional 
abundances by definition add up to 1.0 in the mass isoto- 
pomer envelope analyzed, whatever the value for p or £ The 
absolute value of the enrichments i+l^A^, or EMJ plus deple- 
tions (— A/4^, or DMx), or the degree of perturbation from 
natural abundance values (A/1"), increases as por /increases, 
but EMx is always cancelled out by DM^. 

Conversion of Experimental Data into Biosynthetic 

Parameters 

The sequence followed in converting experimental data 
into biosynthetic parameters is as follows. The ratio (R) of 
Afractional abundances of mass isotopomers in a mixture of 
polymers is used to calculate p and A/1" by using the 
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Table 7. Theoretical abundances of selected mass 



isotopomers in peptide SVVLLLR as a function ofp 



p 














0.00 


0.97702 


0.02293 


0.00005 








0.01 


0.94800 


0.05099 


0.00102 


0.02806 


0.00096 


0.03427 








0.00252 


0.05498 


0.00247 


0.04485 


0.03 


0.89175 


0.10371 


0.00455 


0.08078 


0.00449 


0.05564 


0.04 


0.86451 


0.12840 


0.00709 


0.10547 


0.00703 


0.06667 


0.05 


0.83787 


0.15202 


0.01011 


0.12909 


0.01006 


0.07793 


0.06 


0.81181 


0.17458 


0.01362 


0.15165 


0.01356 


0.08944 


0.07 


0.78633 


0.19609 


0.01758 


0.17317 


0.01753 




0.08 


0.76142 


0.21659 


0.02198 


0.19367 


0.02193 


0.11323 


0.09 


0.73709 


0.23610 


0.02681 


0.21317 


0.02676 




0.10 


0.71332 


0.25463 


0.03205 


0.23170 


0.03200 


0 13811 


0.11 


0.69011 


0.27220 


0.03769 


0.24927 


0.03763 


0.15097 


0.12 


0.66746 


0.28884 


0.04370 


0.26591 


0.04365 




0.13 


0.64535 


0.30457 


0.05008 


0.28165 


0.05003 


0 17762 


0.14 


0.62378 


0.31942 


0.05681 


0.29649 


0.05675 


0 19142 


0.15 


0.60274 


0.33339 


0.06387 


0.31046 


0.06382 


0.20556 


0.16 


0.58223 


0.34652 


0.07126 


0,32359 


0.07120 


0.22004 


0.17 


0.56223 


0.35882 


0.07895 


0.33589 


0.07889 


0.23488 


0.18 


0.54275 


0.37032 


0.08693 


0.34739 


0.08688 


0.25009 


0.19 


0.52376 


0.38104 


0.09520 


0.35811 


0.09514 


0.26569 


0.20 


0.50528 


0.39099 


0.10373 


0.36806 


0.10368 


0.28169 


0,21 


0.48728 


0.40020 


0.11252 


0.37728 


0.11247 


0.29810 


0.22 


0.46975 


0.40869 


0.12155 


0.38577 


0.12150 


0.31495 


0.23 


0.45270 


0.41648 


0.13081 


0.39355 


0.13076 


0.33225 


0.24 


0.43611 


0.42359 


0.14030 


0.40066 


0.14024 


0.35003 


0.25 


0.41998 


0.43003 


0.14998 


0.40711 


0.14993 


0.36829 



Best-fit equation for calculating p.p= -0.0239 + 1.0236(A>i;/A>iJ) - 
1.0414 {^Ayl^AlY + 1.0160(A^6/A/43)3 - 0.8232(A/l|/A/43)'' + 
0.3783A/l6/A/i3)'^. Best-fit equation for calculating as a function 
of p. ^A% = 1.8234(p) - 3.694{pZ) + 2.0276(p3) -i- l.0549{/j4) - 
1,761 HP'S). 



appropriate reference table (e.g., Table 7), AA^ represents the 
A/4;j.(mixture) for 100% newly synthesized molecules at that 
value of p, i.e., the asymptotic A/4^(mixture) that could be 
attained at that value of p. The observed Afractional abun- 
dance of any or all convenient mass isotopomers in the 
mixture [A/l;t(mixture)] is then compared with the maximal or 
asymptotic value, A/1", to generate f. It should be noted that 
both p and Ay4" are directly inferred from R in the reference 
table. That is, in practice one does not interpose an estimate 
of p to calculate A/4^ from experimental data, but A/1^ is 
derived directly from R in the reference table. If an incom- 
plete ion spectrum is monitored (as, in this example, Mq, M3, 
Ms), the correction factor equation is applied for calculating / 
{Eg. A9b). 

Sample Calculations 

A sample of experimental data will now be presented. The 
example described is of an oligopeptide fragment of human 
serum albumin (SWLLLR, or serine-valine-valine-leucine- 
leucine -leucine-arginine), which can be produced by partial 
proteolytic degradation of human serum albumin by treat- 
ment with the enzymes trypsin and chymotrypsin (28). This 
oligopeptide contains three repeating leucine subunits and 
was synthesized in vitro from mixtures of natural abundance 
leucine and [^HsJleucine subunits. The reference table for this 
compound, establishing the relationship between pand Afrac- 
tional abundances, must first be generated. The process is as 
described above. Briefly, the chemical composition of the 
molecule is established (C37N10O9H71). The monoisotopic mass 
of this molecule is 799.6. The molecule is then divided into a 
constant moiety (CsyNioOgHgz) and a variable moiety (*H9). 
The fractional abundance distribution of mass isotopomers 



for the constant moiety is calculated (all elements at natural 
abundance), and the fractional abundance distribution of 
mass isotopomers for each combination of labeled and natural 
abundance subunits in the variable moiety (Hg^Ho, Hg^Hs, 
Ha^Hg, Ho^Hg) is calculated (Table 6/1). Then the fractional 
contribution to the peptide of each, combination of labeled and 
unlabeled subunits is calculated at a given value of p (Table 
6S). For p = 0.10, the likelihood of HgZHo (zero units of 
labeled ^Hj + 3 units of natural abundance labeled H3) = 
0.729, representing the proportion of Mg in the variable 
moiety; the likelihood of Hg^Hs (1 labeled unit + 2 natural 
abundance labeled H3 units) = 0.243, representing the propor- 
tion of M12 in the variable moiety, and so on. Each of these 
coefficients is then multiplied times the mass isotopomer 
fractional abundances in the corresponding combination of 
labeled and natural abundance H3 subunits and incorporated 
into a table of the final abundances of each mass isotopomer 
in the oligopeptide at p = 0.10 (Table QB). The same process 
may be iterated for p = 0.00 to p = 0.20, at stepwise intervals 
of 0.001, for example (Table 7). The natural abundance values 
are then subtracted for each mass isotopomer at each value of 
p to generate Afractional abundances (A/l" or EM^, Table 7). 
The result is a reference table for Afractional abundances of 
mass isotopomers of SWLLLR as a function of p ([^HsJleucine 
proportion in the precursor pool). 

Experimental data from a mixture of isotopically perturbed 
and natural abundance SVVLLLR molecules can then be 
analyzed. The peak areas are converted to fractional abun- 
dances, assuming measurements are made on the complete 
ion spectrum {Mq - Mg). Baseline (natural abundance) experi- 
mental data are subtracted to generate Afractional abun- 
dances. If, for example, the value of AA3 (mixture) = 0.0941 
and AAg (mixture) = 0.0214, then the ratio A/l6(mixture)/ 
A/l3(mixture) was 0.227. The reference table equation (Table 
7) indicates that p = 0.165 and AA3 = 0.2091. Next, f is 
calculated as A/l3(mixture)/A/43 = 0.0941/0.2091, or 45%. In 
a biological experiment, the half-time {tm) of albumin could 
then be calculated (i.e., /= 1 - e*»', and t^z = 0.69/As). 



A ppendix Equations 
Equation A i 

A/! "Sl 

f{X^,%/l^^. . ."-'I,.,,.,) = /^Ab,A-,A-2 . . . x„.,) = ;;r7 n ^? 



where ;^/is the number of atoms of the '/isotope, Nis the total 
number of atoms (A^= S^Jp' -^j). is the number of isotopes of 
/, in this case three (°1, and ^7), and Aj is the fractional 
abundance of '/, where Aj > 0 for i = 0, 1, 2...n-l and 
2"=o Ai= 1. Fis the probability (or statistical frequency) of 
observing the combination of events indicated. Taldng acetic 
acid as an example 
Equation A2a 

/^'3ci2C, 'H4. '6O2) = [pC^C'^O X pOHJ X p('602)] 

= [p(M,C2) X p{.lVk,W,) X p(H.02)] 
= (0.02156)(0.99938)(0.99519) 
= 0.021443 
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Equation A2b 

= [pCMo.CJ X p(M,H4) X p(Mo.02)] 
= (0.97832) (0.00062) (0.995 19) 
= 0.000604 

Equation A2c 

H''C2, 'H4, '^O'^O) = WC,) X PCH4) X pO^O^'O)] 

= [p(Mo,C2) X p(H.H4) X pCM.Oj)] 
= (0.97832) (0.99938) (0.00074) 
= (0.000724) 

where p refers to the probability of that particular element 
group combination occurring, and where p{Mg,Ii) refers to the 
probability of in the element group that follows. 
Equation A2d 

= 0.021443 + 0.000604 + 0.000724 = 0,02277 

A more systematic approach for calculating mass isoto- 
pomer distributions is by disaggregating the molecule into its 
component element groups (e.g., C's, H's, O's), calculating the 
mass isotopomer distributions for each group, and then 
multiplying the appropriate components of each group to- 
gether, systematically, to generate the mass isotopomer distri- 
bution of the molecule 

Equation A3 

U p{.M„.AB) = ^p{,M„.A)Xp{M..m 

n = 0 /=0 " 

A and B can be either element group or molecular distribu- 
tions, and the new, combined distribution is given by the set of 
mass isotopomer fractional abundances of through 
M(„^+„^_ 1). In this way an entire molecular distribution can be 
systematically calculated using Eq. A3 and element group 
mass isotopomer distributions. 

The abundance distribution of labeled and unlabeled sub- 
units that combine to form a polymer is calculated using the 
multinomial expansion 

Equation A4 

Abundance distribution of subunits in a polymer 

= d{z, a, p) = y (^^^ (1 - p)'^-V 

where a is the number of labeled subunits present in the 
variable moiety of the polymer, z is the maximum number of 
subunits that can be labeled in the variable moiety of the 
polymer, and p is the fraction of isotopically labeled subunits 
in the subiinit precursor pool. 

An important mathematical feature of simple mixtures of 
molecular populations is that they combine linearly 

Equation A5a 

Ai in mixture = /(y4, j -I- (1 - /)(A, j 
Equation A5b 



where fis the fraction of polymer population a present (e.g., 
newly sjmthesized polymers), (1 - ^ is the fraction of polymer 
population b present (e.g., preexisting, unlabeled polymers), 
and A^^ and A^ represent the fractional abundances of mass 
isotopomer A^in populations a and b, respectively. 

In a mixture of natural abundance and isotopically per- 
turbed polymers, the fractional abundance of a mass isoto- 
pomer Ax is therefore (the general case in Eqs. A5a and A5b 
above) 

Equation A5c 

A^ (mixture) = {[A,^ + A/l^ + (1 - {){A,)=A,^ + /• kA^ 

where A^^ is the natural fractional abundance of the mass 
isotopomer M„ and is the change in fractional abundance 
of the mass isotopomer Mj, in the isotopically perturbed 
poljoners. We have elsewhere used the term EM^ in place of 
or LAx (13, 14, 16, 17); the terms are interchangeable 
when A/1" is a positive number. It should be noted that A/1" 
and ti.Ax may be positive or negative numbers. 

Equation A5c restates the dependence of A^ on both p (since 
p determines A/4" in the isotopically perturbed polymers) and 
f, in addition to explicitly recognizing the contribution from 
natural abundance molecules {A^^ to the fractional abun- 
dance of the mixture. 

Similarly, the fractional abundance of another mass isoto- 
pomer, Ay, in a mixture is 

Equation A5d 

(mixture) = f(,Ay^+ A/lp + (1 - ti{Ay)=Ay^ + /• A/i; 

where Ay^ is the natural fractional abundance of the mass 
isotopomer My, and A/1" is the change in fractional abundance 
of A^in isotopically perturbed polymers. 

Because of the contributions from natural isotope abun- 
dances (i4j^ and Ay^, the mathematical relationship between 
Ax and Ay in the mixture (the internal pattern or distribution) 
is not uniquely determined by p, which determines the mass 
isotopomer pattern in the isotopically perturbed molecules, 
but is also determined by / 

Equation A6 

Ax (mixture)M^ (mixture) = [A^^ + f- ^AyJ/{Ay^ + /• AAp 

Accordingly, there are two unknowns that influence the 
internal relationships among mass isotopomers in a mixture 
of labeled and natural abundance polymers. More than one 
equally valid mathematical approach has been proposed to 
solve for p and f from experimental data in the presence of 
natural abundance contributions (6, 12, 14-17, 20, 21, 32). 
Our approach has been to linearize the relationship between 
p and /l;i.(mixture) by expressing the results as Afractional 
abundances (excesses or depletions compared with natural 
fractional abundance values). When Afractional abundances 
are used, the internal relationship among mass isotopomers 
in the mixture [A/l^(mixture)/A/4y(mixture)], when the changes 
in abundance are in the positive direction, becomes constant 
for p and independent of f, so that the initial reference table 
can be used regardless of the particular mixture present 

Equation A7a 

AAx (mixture) = -I- /• A/1^] - [A,^] = f- AA'^ 
Equation A7b 



A2 in mixture = fiA^) + (1 - UAi) 



LAy (mixture) = \Ay^ -(- /• A/ip - [^J = /• A^; 
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Equation A 7c 

(mixture)/A/ly (mixture) = /• ^A*lf■ A.4p = A^"/A.4^ 

Because Ay4^ and A/1" are defined as the change in frac- 
tional abundance in the newly synthesized or isotopically 
perturbed polymers only, Eq. A7c is independent of the 
dilution by unlabeled polymers, and thus only depends on the 
value of p that was present. The value of p is then calculated 
from the best-fit equation of jD vs. ^A'^/AA'^ in the appropriate 
reference table. Using the example of the peptide SVVLLLR, 
HAAe/AAs = 0.19142 (Table 7) 

Equation A8 

p(SWLLLR) = -0.0239 + 1.0236(0.19142) 

- 1.0414(0.19142)2 -I- 1.0160(0.19142)3 _ 0.8232(0.19142)^ 

-I- 0.3783(0.19142)= = 0.1400 

The value of /is then calculated using AA", at the value of p 
calculated from the AAJAAy ratio present 
Equation A9a 

f= LA^ (mixture)/AA" 

where A/4" represents the asymptotic value of Ai4;,(mixture) 
(i.e., if 100% of molecules were labeled at the p calculated), 
corrected for incomplete ion spectrum sampling, if necessary 
(see further discussion). 

When all of the ions in the mass isotopomer spectrum are 
not monitored, there is no effect on calculation of />, but an 
effect on /may be observed, unless the calculation algorithm 
is corrected for the proportion of ions sampled from the 
labeled and unlabeled populations in the mixture of polymers 
(see text of article). Because of unequal molar contributions to 
the ion spectrum sampled, a mole of natural abundance 
molecules will be overrepresented in the envelope sampled 
relative to a mole of enriched molecules, and the assumption 
of linearity for mixtures (Eq. A7a) no longer applies. It is 
therefore necessary to weight the contributions from different 
populations. A simple algebraic correction can be formulated 

If 

R{£) = ratio or proportion of ions in the total envelope 

monitored, at baseline 
R(£) = ratio or proportion of ions in the total envelope 

monitored, for enriched molecules 
M(meas) = measured Mi fractional abundances within the 

spectrum sampled 
M\ (B) = True (theoretical) proportion of Mj ions in the 

total envelope, at baseline 
M(£) = True (theoretical) proportion of Mi ions in the 
total envelope, for enriched molecules, then: 
Equation A9b 



[R(B) X M(^eas)l - [Mi(B)] 
[M,(E) - Mi,B)] + [(R(B) - R(E)) X M,(„,^,] 

The derivation is as follows 

[/. M„£,] + [(1 - 4 X Mi(B,] 



ARe - Rb) + Rs 

- M,(B)1 



when R(jg = R(bi = 1.0 (i.e., all ions 
populations), the equation reduces to 



+ [(Rb - Ri X A^I(„,eas)] 

monitored, for both 



- Mi(^ AAi (mixture) 



So the equation reduces to the standard form (Eq. A9a), as 
it should, when the complete ion spectrum is monitored. 

Address for correspondence and reprint requests; M. K. Heller- 
stein, Dept. of Nutritional Sciences, Univ. of California, Berkeley, CA 
94720 (E-mail: march@nature.berkeley.edu). 
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